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ABSTRACT 

Non-coherent  detection  of  Gaussian  targets  (Swerling  II  targets)  in  the  K-distributed  clutter 
environment  is  investigated.  The  optimal  detector  is  derived  based  on  the  Neyman-Pearson 
principle.  It  is  shown  to  be  the  well-known  square-law  detector.  Amplitude  detector,  log 
detector,  and  the  like  are  not  optimal,  and  result  in  some  detection  loss.  Temporally  correlated 
clutter  provides  a  target  gain,  and  improves  detection.  The  higher  the  temporal  correlation, 
the  higher  the  target  gain.  Spatially  correlated  non-Gaussian  clutter  can  also  provide  a  CFAR 
gain.  The  autoregressive  technique  is  used  to  optimally  estimate  the  texture  of  the  clutter. 
That  in  turn  significanly  improves  the  detection  compared  to  the  traditional  cell-averaging 
processing. 
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Optimal  Detection  in  the  K-Distributed  Clutter 
Environment  --  Non-Coherent  Radar  Processing 


Executive  Summary 

In  the  maritime  environment,  radar  detection  unavoidably  needs  to  deal  with 
undesired  signals,  primarily  sea  clutter  (echoes  from  the  sea  surface).  How  to  detect 
target  signals,  especially  relatively  weak  ones,  against  the  clutter  is  challenging. 

Optimal  detectors,  based  on  the  Neyman-Pearson  principle,  that  maximise  the 
probability  of  detection  for  a  given  false-alarm  rate  are  the  most  interesting  detectors  to 
the  radar  community.  Their  derivation  depends  on  the  statistical  models  of  both  the 
clutter  and  the  target.  In  a  recent  paper,  'Optimal  coherent  radar  detection  in  a  K- 
distributed  environment'  we  have  discussed  the  problem  of  optimal  coherent 
detection.  This  report  focuses  on  the  non-coherent  detection  of  Gaussian  targets 
(Swerling  II  targets)  in  the  compound  K-distributed  clutter  environment. 

This  report  makes  the  following  three  contributions. 

First  the  optimal  detector  for  multi-look  non-coherent  detection  of  Gaussian  targets  in 
the  compound  K-distributed  clutter  is  derived.  The  optimal  detector  derived  is  shown 
to  be  the  well-known  square-law  detector.  This  is  because  the  clutter  undergoes  a 
Gaussian  random  process  during  the  multi-look  processing  period  (i.e.,  the  multi-pulse 
processing  period),  as  the  slowly-varying  component  of  the  compound  clutter  remains 
unchanged  during  the  period,  according  to  the  assumption.  Although  the  derived 
optimal  detector  is  not  new,  the  derivation  itself  has  a  guiding  meaning.  As  the 
detector  has  been  rigorously  derived  for  the  first  time  using  the  Neyman-Pearson 
principle,  it  means  that  no  other  detectors  exist  which  would  perform  better  for  the 
given  conditions.  Other  detectors,  such  as  the  multi-look  amplitude  detector,  the  multi¬ 
look  log  detector,  and  the  like  are  not  optimal  and  inherently  result  in  some  detection 
loss. 

Secondly  we  have  shown  that  for  temporally  correlated  clutter,  the  use  of  a  multi-look 
whitening  process  provides  a  target  gain  and  improves  the  detection.  The  higher  the 
correlation,  the  larger  the  target  gain.  The  target  gain  comes  from  the  difference 
between  the  spectrum  of  the  correlated  clutter  and  the  spectrum  of  uncorrelated  target 
signals  (a  non-uniform  spectrum  against  a  uniform  spectrum),  providing  a  second 
characteristic  (in  addition  to  the  intensity)  for  discriminating  the  target  from  the  clutter. 
On  the  other  hand,  if  both  the  clutter  and  the  target  signals  are  individually 
uncorrelated  (the  cross-correlation  between  the  two  is  always  zero),  each  of  them  has  a 
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uniform  spectrum,  and  there  is  only  one  characteristic  (intensity)  that  can  be  employed 
in  the  detection.  Therefore  the  event  of  the  uncorrelated  Gaussian  targets  embedded  in 
the  uncorrelated  compound  clutter  represents  the  worst  scenario  in  terms  of  detection. 
If  the  correlation  of  the  target  signals  is  the  same  as  the  correlation  of  the  clutter,  a 
treatment  of  de-correlation  leads  to  the  same  processing  for  the  uncorrelated  target  in 
the  uncorrelated  clutter. 

Lastly,  spatially  correlated  non-Gaussian  clutter  may  be  able  to  provide  some  constant 
false-alarm  rate  (CFAR)  gains.  The  CFAR  gain  is  dependent  on  the  estimate  of  the  local 
mean.  For  this  analysis  we  have  examined  the  use  of  the  linear  autoregressive 
technique  and  derived  the  optimal  weights  for  estimating  the  local  mean  of  clutter.  The 
autoregressive  estimation  is  optimal  under  the  linear  assumption  and  better  than  the 
traditional  cell-averaging  estimation.  The  optimal  estimation  (under  the  linear 
assumption)  of  the  clutter  texture  has  in  turn  resulted  in  a  further  significant  detection 
improvement  (a  few  dB)  for  highly  spatially  correlated  K-distributed  clutter  compared 
to  the  traditional  cell-averaging  estimation. 

This  work  was  carried  out  in  support  of  the  ADF's  Air  7000  Project. 
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1.  Introduction 


In  a  journal  paper,  'Optimal  coherent  radar  detection  in  a  K-distributed  clutter 
environment'  (Dong  2012),  we  have  proposed  an  optimal  detector  for  coherent  detection 
against  K-distributed  clutter.  This  report  develops  optimal  and  near  optimal  detectors  for 
non-coherent  detection. 

In  the  past,  radar  systems  had  relatively  low  resolution  capabilities,  and  the  Gaussian 
clutter  assumption  was  a  valid  model  for  the  clutter.  With  the  advances  in  radar 
technology,  the  resolution  capabilities  of  radar  have  been  greatly  improved  in  recent  years. 
Clutter  collected  using  high  resolution  radar  systems  exhibits  non-Gaussian  behaviour. 
The  associated  problem  of  optimal  detection  in  a  non-Gaussian  background  remains  to  be 
solved.  The  compound  K-distribution,  for  instance,  is  a  non-Gaussian  distribution 
commonly  used  to  model  radar  sea  clutter;  one  of  the  milestone  findings  in  radar  clutter 
analysis  and  research  in  resent  years  (Ward  et  al.  2006). 

Researchers  have  been  studying  the  optimal  coherent  detection  of  radar  targets  embedded 
in  compound-Gaussian  clutter  for  many  years  (Sangston  et  al.  2010;  Sangston  and  Gerlach 
1994;  Sangston  et  al.  1999;  Gini  et  al.  1999;  Farina  et  al.  1997;  Gini  et  al.  1998).  While  the 
matched  filter  is  optimal  for  Gaussian  clutter,  the  paper  (Dong  2012)  shows  that  the 
proposed  optimal  detector  for  K-distributed  clutter  significantly  improves  the  detection 
compared  to  the  matched-filter  if  clutter  is  highly  spiky  (i.e.,  K-distributed  clutter  with  a 
small  shape  parameter).  Therefore,  there  seems  a  need  to  investigate  the  optimal  non¬ 
coherent  detection  against  the  K-distributed  clutter. 

Non-coherent  detection  against  K-distributed  clutter  has  received  equal  attention. 
Armstrong  and  Griffiths  (1991)  studied  detection  of  fluctuating  targets  in  spatially 
correlated  clutter,  but  their  study  focused  on  the  performances  of  cell-averaging  (CA),  cell 
averaging  greatest  of  (CAGO),  and  ordering  statistic  (OS)  constant  false-alarm  rate  (CFAR) 
processors  and  did  not  discuss  the  issue  of  optimal  detection.  Watts,  Ward  and  Tough 
studied  CFAR  loss  and  CFAR  gains  associated  with  the  K-distributed  clutter  (Watts  1996; 
Watts  et  al.  2007;  Watts  1987).  It  has  been  found  that  for  the  compound  Gaussian  clutter, 
such  as  K-distributed  clutter,  CA-CFAR  can  provide  a  CFAR  gain,  provided  that  the 
texture  of  clutter  (underlying  mean)  is  correlated  (Watts  1985).  The  higher  the  correlation, 
the  larger  the  CFAR  gain.  Watts  (1985)  thus  proposed  a  concept  of  'ideal  CFAR'  which 
means  that  if  the  exact  mean  of  the  clutter  for  the  cell  under  test  (CUT)  is  known  or 
estimated  by  other  means,  the  best  performance  can  be  achieved.  Based  on  this  concept, 
Buccuarelli  et  al  (Bucciarelli  et  al.  1996)  proposed  to  use  the  maximum  a  posteriori 
estimation  in  the  logarithm  domain  (LMAP)  to  estimate  the  local  mean  for  the  CUT.  They 
found  that  the  performance  of  LMAP-CFAR  outperforms  CA-CFAR  especially  when  the 
correlation  of  the  texture  is  high. 

A  pulsed  Doppler  radar,  equipped  with  a  single  transmitter  and  a  single  receiver,  usually 
collects  two-dimensional  data,  one  dimension  is  time  (separated  by  pulses)  and  the  other 
is  range  (separated  by  range  bins).  This  kind  of  data  is  often  referred  to  as  multi-look  data. 
The  associated  detection  problem  is  how  to  process  this  two-dimensional  data  to  achieve 
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the  detection  goal.  Apparently  the  above  studies  of  non-coherent  detection  only  consider 
the  processing  in  the  range  domain. 

Conte  et  al  (1999)  proposed  a  generalised  likelihood  ratio  test  (GLRT)  model  for  the 
optimal  incoherent  detection  of  Swerling  II  targets  in  the  compound  K-distributed  clutter. 
However  the  model  assumes  that  for  the  1 1  ]  hypothesis  only  the  first  portion  of  an 
incoherent  pulse  train  contains  the  clutter  and  target  signal,  and  the  second  portion  the 
pulse  train  contains  clutter  only.  This  condition  does  not  seem  to  be  robust  to  an  unknown 
target  location. 

In  this  report  we  first  consider  the  processing  in  the  time  domain.  Different  assumptions 
for  target  signal  and  clutter  lead  to  different  designs  of  detection  scheme  in  order  to 
achieve  the  optimal  and  near-optimal  performance.  This  report  considers  Swerling  II 
model  as  the  target  model,  i.e.,  target  RCS  is  independent  from  pulse  to  pulse  and  varies 
with  Gaussian.  For  such  a  radar  target,  coherent  detection  is  not  appropriate,  because 
target  signal's  phase  varies  randomly,  and  both  target  signals  and  clutter  have  a  uniform 
spectrum  in  the  frequency  domain.  Often  we  need  to  detect  the  intensity  of  the  target 
signals.  The  intensity  is  enhanced  by  multiple  look  non-coherent  integration.  For  Gaussian 
clutter,  the  multi-look  intensity  averaging  processing  is  optimal.  However  whether  it  is 
still  optimal  for  non-Gaussian  clutter  is  remains  to  be  answered.  This  report  tries  to  derive 
optimal  and/or  near  optimal  detector  from  the  Neyman-Pearson  principle  (Kay  1998, 
page  174). 

The  unwanted  signals  that  need  to  be  considered  in  the  maritime  environment  include 
echoes  from  the  sea  surface  plus  thermal  noise  of  the  radar  receiver.  In  this  report,  the 
unwanted  signals  are  assumed  to  be  represented  by  a  compound  K-distribution  (this  is 
justified  in  Section  2).  We  derive  optimal  and  near-optimal  detectors  based  on  the 
Neyman-Pearson  principle  in  Section  3.  The  performance  of  optimal  and  near  optimal 
detectors  for  temporally  correlated  clutter  is  also  analysed. 

Secondly  we  study  the  processing  in  the  range  domain.  Specifically  we  consider  in 
Section  4  the  non-coherent  detection  for  the  spatially  correlated  clutter.  We  use  the 
autoregressive  (AR)  technique  to  optimally  estimate  the  texture  of  clutter.  That  in  turn 
significantly  improves  the  detection  compared  to  the  traditional  CA  processing. 


2.  Justification  of  K-distributed  Clutter 

In  a  maritime  radar  surveillance  environment,  sea  clutter  has  been  verified,  through 
numerous  trials,  to  fit  with  the  compound  K-distribution  for  most  conditions  (Ward  et  al. 
2006;  Crisp  et  al.  2006;  Dong  and  Merrett  2010;  Greco  and  Gini  2007;  Farina  et  al.  1997). 
Horizontally  polarised  higher  resolution  sea  clutter  may  even  have  a  heavier  tail  in  its 
probability  density  function  (pdf)  and  fits  better  with  other  distributions,  such  as  KA 
(Watts  et  al.  2005),  KK  (Dong  and  Haywood  2007)  and  Pareto  (Farshchian  and  Posner  2010) 
distributions.  This  report  assumes  sea  clutter  to  be  K-distributed.  The  compound 
K-distribution  is  composed  of  two  components,  a  fast-varying  component,  referring  to 
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speckle  and  a  slowly- varying  component,  referring  to  the  underlying  mean  or  texture.  The 
fast-varying  component  is  a  zero  mean,  unit  variance  complex  Gaussian  process 
modulated  by  the  slowly-varying  component,  whose  intensity  is  gamma  distributed 
(Ward  ef  al.  2006). 

Thermal  noise,  according  to  its  nature,  is  often  modelled  by  a  Gaussian  random  process. 
The  distribution  of  the  sum  of  the  K-distributed  clutter  and  Gaussian  thermal  noise 
unfortunately  does  not  have  closed-form.  To  overcome  this.  Watts  (1987)  has 
approximated  the  combined  distribution  as  a  new  K-distribution,  by  equating  their 
intensities  and  variances.  In  other  words,  the  combined  distribution  of  a  K-distribution 
with  parameters  of  juc  (mean  intensity)  and  vc  (shape  parameter)  and  a  Gaussian  thermal 

noise  with  parameter  <j~  (variance  or  mean  intensity)  is  approximated  as  a  new  K~ 
distribution  with  parameters  of  fi  and  v .  The  new  parameters  /u  and  v  are  determined 
by  equating  the  mean  and  the  variance  of  the  two  distributions  (Watts  1987).  Through 
numerical  simulation,  we  found  that  this  approximation  provides  very  good  agreement 
between  the  theoretical  and  data  distributions  for  a  wide  range  of  Gaussian  signals  from  a 
very  weak  thermal  noise  (very  small  signal-to-clutter  ratio  (SCR))  to  very  strong  Gaussian 
target  signals  (very  large  SCR)  (see  Appendix  B  for  details). 

Sea  clutter  data,  whose  true  distribution  is  unknown,  when  received  by  a  radar  system 
includes  the  thermal  noise  of  the  radar.  Consequently  the  estimated  mean  and  shape 
parameters  will  automatically  be  for  the  combined  distribution.  Therefore,  the  undesired 
signals  considered  in  this  report  only  have  a  single  component,  the  combined  compound 
K-distributed  clutter. 


3.  Optimal  Detector  for  Non-Coherent  Detection 

3.1  Uncorrelated  Clutter 

In  order  to  employ  the  Neyman-Pearson  principle  to  indentify  the  optimal  detector  we 
need  to  generate  a  statistical  representation  of  the  target  embedded  in  the  clutter 
distribution.  The  mathematical  analysis  for  the  case  under  consideration,  a  Gaussian  target 
embedded  in  K-distributed  clutter  plus  noise  is  provided  in  Appendix  B. 

Having  approximated  the  distribution  of  Gaussian  target  embedded  in  K-distributed 
clutter,  we  are  now  in  a  position  to  discuss  the  associated  optimal  detector  for  non¬ 
coherent  detection.  The  multi-look  case  means  that  multiple  measures  jc[«],  n  =  1 ,  •  •  ■ ,  N , 
with  respect  to  pulse,  for  each  range  bin  are  available.  We  want  to  derive  a  detection 
scheme  which  is  optimal,  i.e.,  the  probability  of  detection  is  maximum  for  a  given  false- 
alar  m  rate. 

The  Neyman-Pearson  principle  (Kay  1998,  page  174)  states  that  to  maximise  a  probability 
detection  Pd  for  a  given  probability  of  false-alarm,  the  detection  threshold  is  based  on  the 
likelihood  test  ratio  of. 
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Pte  Ho) 


(1) 


0>  H  o 


where  H0  and  //,  are  null  and  alternative  hypotheses  of  absence  and  presence  of  target 
signal.  We  assume  that, 

(a)  the  clutter  parameters,  namely,  the  mean  and  the  shape  parameter,  are  known,  or 
can  be  estimated  from  the  secondary  data; 

(b)  during  the  multi-look  sampling  period  the  slow-component  of  the  K-distributed 
samples  remains  unchanged;  and 

(c)  multi-look  samples  are  independent  (the  correlated  case  is  discussed  later). 

Since  the  distribution  under  //,  is  approximated  by  another  K-distribution  (see  Appendix 
B),  the  marginal  probability  density  function  (pdf)  is  (Sangston  et  al.  2010), 


<  =  o,i 


o  7i  r. 


(2) 


where  p( r)  is  the  gamma  distribution  of  the  slowly-varying  underlying  component: 


=  exp(-  br) 

r(v) 


(3) 


where  v  is  the  shape  parameter,  b  =  — ,  ju  =  E{\  x\2}  is  the  intensity  mean.  The  integral 

M 

gives. 


p(x;Hi)  = 


x"nv,) 


Vj-N  Vj+N 
2 


A  2  b.  2  K  „  2  Jb.A  i  =  0,l 


(4) 


H  N 

where  N  is  the  number  of  multi-looks  and  A  =  x  x  =  X I  x[n\  I"  •  Finally  we  have, 

n= 1 


In  A(x)  =  lnr(v0)  -lnr(Vj) 


2  2 

In  KVi_n  (2  VM  )-  ^  KVq_n  (2  7m  ) 


V-^\nA  +  V-^\nb,-V-^\nb0 


(5) 
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The  likelihood  ratio  test  (LRT)  of  (5)  is  a  non-linear  function  of  measure  A  and  the  design 
of  the  optimal  detector  based  on  (5)  is  difficult  as  the  pdf  of  In  a(\)  with  respect  to 
measure  A  is  unknown.  Next  we  examine  how  the  LTR  of  (5)  behaves  and  consequently 
construct  a  near-optimal  detector. 

Since  the  marginal  pdf  p(x;Hi)  is  only  a  function  of  measure  A  ,  if  we  replace  p(x;  //  ) 
by  p(z;Hj)  ,  we  find  that  the  corresponding  LRT  of  (5)  remain  unchanged,  where 
p{z\H.)  is  the  pdf  of  multi-look  K-distributed  data  given  by  (B6)  in  Appendix  B,  and 
z  =  N~l  X,)=i  I  x[«] !  -AIN.  The  advantage  of  using  this  replacement  is  that  the  pdf  of 
p(z;  IF  )  can  be  plotted  against  z  .  In  another  word,  the  LRT  is  replaced  by. 


A(z) 


P&  gl)  I1  y 

P(N  H0)  h0 


(6) 


One  may  verify  (6)  and  (1)  are  identical,  i.e..  In  A(x)  —  In  A(z) . 

Figure  1  shows  the  LTR  In  A(z)  varies  against  the  measure  z  given  parameters  of  v  =  1.2, 
p  — 1.0,  N  =  10  and  SCR—  6 dB .  Initially,  it  seems  difficult  to  understand  the  curve.  For 
instance,  if  one  selects  In  A (z)  =  0  as  a  threshold  (the  selection  of  the  threshold  only 
depends  on  Pfa ),  then  the  optimal  detector  only  declares  presence  of  target  when  the 
measure  z  falls  in  a  certain  region  (between  3.73dB  and  14.81dB).  It  is  easy  to  understand 
that  the  optimal  detector  rejects  //,  when  measure  z  is  low  (z  <  3.73 dB ).  However,  why 
does  the  optimal  detector  also  rejects  Hx  when  measure  is  high  (z  >  14.81  dB )?  The  next 
paragraph  explains. 


Figure  1:  LTR  In  A(z)  against  measure  z  =  N  1  X„'=i  |  *[«]  f 
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Shown  in  Figure  2  are  pdfs  of  multi-look  clutter  and  multi-look  target  plus  clutter, 
respectively,  for  the  given  parameters.  The  pdf's  are  shown  on  linear  scale  and  log  scale, 
respectively.  Shown  in  the  figure  is  also  the  LTR.  Now  the  interpretation  of  LTR  curve 
becomes  clear.  When  z  value  is  low,  p(z;H0 )  is  higher  than  p(z;  //, ),  and  the  difference 
gradually  decreases  with  increasing  in  z  ,  and  reaches  a  point  Tx  ,  where 
p(T{;H0)  =  p(T] ;  //, ) .  On  the  log  scale,  LTR  monotonically  increases  and  reaches  0  at  ZJ . 
After  T{ ,  p(z;  If )  becomes  higher  than  p(z\  //0),  so  LTR  continues  its  monotonic  increase 
and  reaches  its  peak  where  the  ratio  of  the  two  is  maximum.  After  that  point,  though 
p(z;Hl)  is  still  higher  than  p(z;  //„),  the  ratio  decreases.  The  LTR  curve  starts  dropping 
from  its  peak  till  reaching  0  at  T2  where  again  p(T2;H0)  =  p(T2\Hf).  After  T2  , 
p(z;H0) again  becomes  higher  than  p(z;  //, ),  as  the  distribution  of  p(z;  //,,)  has  a  longer 
and  heavier  tail,  according  to  the  given  conditions.  Therefore,  the  curve  LTR 
monotonically  decreases  after  it  surpasses  its  peak. 


Measure  z,  dB 


Figure  2:  (top  and  middle)  pdfs  of  clutter  and  target  plus  clutter  on  linear  scale  and  log  scale, 
respectively;  (bottom)  LTR.  The  region  satisfi/ing  LTR  greater  than  In  A(z)  =  Ois  also 
shown  in  all  three  plots. 
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Supposing  we  choose  a  threshold.  In  X (z)  >  0 ,  to  accept  //,  as  shown.  This  is  equivalent 
to  accepting  Hl  when  p(zl;Hl)  >  p(z;H0).  Therefore,  the  region  [7j,  T2]  that  the  optimal 
detector  accepts  Hl  starts  at  the  point  where  the  pdf  of  clutter  plus  target  is  higher  than 
the  pdf  of  clutter  and  ends  at  the  point  where  the  pdf  of  clutter  becomes  higher  than  the 
pdf  of  target  plus  clutter.  It  can  be  seen  that  the  optimal  detector  does  what  exactly  is 
required,  because  the  condition  p(zl',Hl )  >  p{z\  H0)  does  not  hold  for  A  <  z  <o o  . 

If  we  relax  the  optimal  condition  from  Tj  <  T  <  T1  to  T  >r,  to  form  a  sub-optimal  detector 
for  accepting  Hl ,  we  need  to  study  the  consequences.  Replacing  7j  <  T  <  T2  by  T  >  Tj  will 
theoretically  increase  both  the  false-alarm  rate  and  the  probability  of  detection,  because, 

Pfa  =  \p{z’,H0)dz  <  ]p(z;H0)dz  (7) 

Ji  TX 

Pd  =  j p(z;  Hx  )dz  <  ] p(z ;  Hx  )dz  (8) 

n  tx 

Our  concern  is  only  the  increase  of  the  false-alarm  rate.  For  the  optimal  detector  accepting 
//,  for  Tx  <  T  <  L, ,  the  corresponding  false-alarm  rate  is, 

pfa  =  jp{z;H0)dz  (9) 

n 

Alternatively,  the  sub-optimal  detector  accepting  //  for  T>TX  accompanies  a  false-alarm 
rate  of, 

]p(z;H0)dz  =  j  p(z;  H  0)dz  +  ]p(z;H0)dz  =Pfa  +  A Pfa  (10) 

Ji  Ji  Ji 

In  order  to  have  a  reasonable  value  of  Pd  for  the  optimal  detector,  the  interval 
P  -  7J  »  0,  as  Pd  =  Qp(z;  II  x)dz .  Because  residuals  of  the  pdf  for  z  >  T2  is  very  small, 
we  generally  have  A Pfa  «  Pfa  ■  For  instance,  for  the  above  example,  we  found, 

Pfa  =0.08754  and  APfa  =  1. 14  x  10_y  ,  and  correspondingly,  Pd  =  0.9552  and 
A Pd  =6.24xl0_l°.  Therefore,  relaxing  the  condition  from  Tx  <  T  <  T2  to  T  >  Tx  does  not, 
in  practical  terms,  change  either  Pfa  or  Pd  for  this  example. 

The  value  of  A Pfa  will  surely  depends  on  the  parameters.  However,  for  a  common  radar 
detection  problem,  it  often  requires  a  small  Pfa  (usually  1 0~"  or  smaller),  and  a  reasonable 
Pd  (usually  0.5  or  higher).  In  order  to  achieve  \^p{z\Hx)dz  =  P,  >0  .5 ,  the  interval  T2  -  Tj 
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must  not  be  too  small.  On  the  other  hand,  the  value  of  p(z;H0)  in  this  region  has  to  be 
small  to  achieve  |r2 p(z;H0)dz  =  Pfa .  As  shown  in  the  above  example,  for  such  conditions, 
«  Pfa  holds.  For  instance,  the  cumulative  distribution  function  (cdf)  of  clutter  for  the 
above  example  is  shown  in  Figure  3.  According  to  this  cdf,  to  achieve  a  false-alarm  rate  of 
10“6,  the  sub-optimal  detector  requires  z  >  12.53d/!,  Pr(z  >  \2.53dB;  H0 )  =  I  0-6  .  On  the 

other  hand,  the  optimal  detector  might  require  12.53(75  <  z  <  \5.52dB  to  have  a  3dB 
interval  to  achieve  a  reasonable  probability  of  detection  for  a  scenario.  However,  the 
relaxation  from  the  optimal  condition  of  12.53(75  <  z  <  15.52(75  to  the  suboptimal 
condition  of  12.53 dB  <  z  <  <x>  have,  in  a  practical  term,  the  same  false-alarm  rate  of  10_(1 
because  of  Pr(z  >  15.52 dB',H0)  <  10“9 . 


Clutter  intensity  (dB) 

Figure  3:  Cdf  of  multi-look  K-distributed  clutter  shown  as  log10(l -oT/")  (v  =  1.2  ,  p-l.O, 
A  =  10). 


Therefore,  the  sub-optimal  detector  using  z  >Tt  for  detection  performs,  for  practical  cases, 
identical  to  the  optimal  detector  using  T{<  z  <T2  provided  that  the  condition  of 
7)  —  7j  »  0  is  satisfied. 

Since  the  sub-optimal  detector  only  requires  z  >  7j ,  it  means  that  the  associated  LTR  curve 
can  be  assumed  to  be  monotonic  for  0  <  z  <  oo .  One  of  such  LTR  curve  is, 

A(z)  =  z  (11) 

The  calculation  of  the  sub-optimal  detector  will  become  tractable,  as  the  pdf's  of  p(z;  //„) 
and  p(z;Hf)  are  known.  As  explained  and  for  simplicity,  the  sub-optimal  detector  using 
the  LTR  of  (11)  will  simply  be  referred  to  as  the  optimal  detector. 
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The  optimal  detector  of  (11)  applies  to  any  shape  parameter  values.  When  the  shape 
parameter  approaches  infinity,  it  becomes  the  case  of  a  Gaussian  target  embedded  in 
Gaussian  clutter.  It  is  well-known  that  the  LTR  of  (11)  is  indeed  the  optimal  for  such  a 
case.  Since  it  is  assumed  that  for  the  K-distributed  clutter  the  underlying  mean  remains 
constant  during  the  multi-look  processing  period,  the  variation  of  x[«],  n=\,  is 

therefore  in  fact  a  Gaussian  random  process  for  the  processing  period.  This  explains  that 
the  detector  determined  by  the  LTR  of  (11)  is  also  optimal  (in  practical  sense)  for  the  K- 
distributed  clutter  whose  slowly-varying  component  (underlying  mean)  is  unknown  but 
remains  unchanged  during  the  multi-look  processing  period.  Since  we  also  assumed  that 
the  underlying  mean  is  spatially  (from  range  bin  to  range  bin)  uncorrelated,  and  the  best 
estimate  of  it  therefore  is  the  global  mean.  The  case  of  spatially  correlated  clutter  will  be 
discussed  in  Section  4. 

Although  the  above  derivation  results  in  an  identical  optimal  detector  for  Gaussian  clutter 
and  the  compound  K-distributed  clutter.  It  has  a  guiding  meaning  in  research.  While  the 
performance  of  the  multi-look  detector  against  K-distributed  clutter  has  been  discussed  by 
many  researchers,  here  we  first  mathematically  prove  that  it  is  also  practically  optimal  for 
the  compound  K-distributed  clutter  using  the  Neyman-Pearson  principle.  It  shows  that 
detection  must  use  the  intensity  data  in  order  to  obtain  the  optimal  detection.  Detection 
schemes  using  other  data  formats,  such  as  amplitude,  logarithm  transform  and  so  on  are 
not  optimal,  and  will  result  in  some  detection  loss  for  the  same  false-alarm  rate. 

Receiver  operating  characteristic  (ROC)  curves1  for  the  optimal  detector  determined  by 
A (z)  =  z  are  shown  in  Figure  4  and  Figure  5  for  various  given  parameters.  It  can  be  seen 
that  for  the  same  probability  of  detection,  the  required  SCR  increases  substantially  with 
the  increasing  in  spikiness  of  sea  clutter,  compared  to  the  exponential  case.  For  instance, 
for  the  case  of  Pfa  =10-6,  Pd  =0.5  and  N  =  10,  it  requires  a  SCR=12.4dB  for  v  =  1.2 

compared  to  a  SCR=3.8dB  for  v  =  oo  (exponential  case).  It  is  understood  that  because  of 
the  spikiness,  the  threshold  has  to  increase  substantially  in  order  to  maintain  the  same 
false-alarm  rate.  The  increase  in  the  threshold  in  turn  substantially  decreases  the 
probability  of  detection.  One  way  to  increase  the  probability  of  detection  is  to  increase  the 
number  of  multi-looks  in  the  non-coherent  processing  if  possible. 


1  Conventionally,  an  ROC  curve  plots  Pd  against  Pfd  for  a  given  SCR.  However,  a  plot  of  Pd 
against  SCR  for  a  given  Pfd  is  also  referred  to  as  an  ROC  curve  recently.  This  report  adopts  the 
latter  definition. 
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Figure  4:  ROC  curves  of  the  optimal  detector  determined  by  A(z)  =  z  with  Pfa 
N  =  4. 


Figure  5:  ROC  curves  of  the  optimal  detector  determined  by  /L(z)  =  z  with  Pfa 
N  =  10. 
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Some  ROC  curves  shown  in  Figure  5  were  confirmed  using  Monte  Carlo  simulation.  In  the 
simulation,  first  a  dataset  of  K-distributed  clutter  and  a  dataset  of  Gaussian  distributed 
target  signals,  both  in  the  complex  domain,  were  generated  according  to  the  given 
parameters.  The  two  datasets  were  combined  together  in  the  complex  domain.  After 
taking  the  square  of  data's  amplitude,  multi-look  average  processing  was  followed.  The 
number  of  measures  exceeding  the  corresponding  threshold  was  counted  and  the 
associated  Pd  was  calculated.  The  detection  performance  calculated  by  the  Monte  Carlo 
simulation  is  shown  in  Figure  6,  together  with  the  correspondingly  theoretical  ROC  curves 
for  comparison.  It  can  be  seen  that  they  are  consistent. 


SCR,  dB 

Figure  6:  Confirmation  of  detection  performance  using  Monte  Carlo  simulation. 


The  performance  of  other  detectors  was  investigated.  In  particular,  we  take  a  look  at  two 
other  detectors,  namely,  the  amplitude  detector  using  the  measure  of  AO'X^il  x[n]  \,  and 

the  log  detector  using  N~l  X,)  i  10  log10  |  x[n]  f  for  detection,  respectively.  While  there  are 
analytical  forms  of  pdf  and  cdf  for  single-look  K-distributed  data  when  the  data  is 
measured  in  amplitude  or  in  log  domain  (a  simple  pdf  transform  governed  by 
p(y(x))  =  p(x)\dx  /  dy\),  the  pdf  and  cdf  of  multi-look  K-distributed  data  do  not  have 

simple  analytical  forms  when  multi-look  average  processing  is  performed  for  the 
amplitude  or  in  the  log  domain.  Therefore  the  performance  of  the  amplitude  detector  and 
log  detector  for  the  Gaussian  target  embedded  in  the  K-distributed  clutter  is  numerically 
calculated  using  the  Monte  Carlo  simulation  (the  accuracy  of  the  Monte  Carlo  simulation 
has  been  demonstrated  in  Figure  6).  Shown  in  Figure  7  are  ROC  curves  of  amplitude 
detector  and  log  detector  in  comparison  with  the  optimal  detector,  i.e.,  the  intensity 
detector.  It  can  be  seen  that  neither  the  amplitude  nor  log  detectors  perform  as  well  as  the 
optimal  detector.  Parameters  used  in  the  simulation  are  given  in  the  figure.  Other 
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parameters  were  also  tested,  and  the  associated  ROC  curves  have  a  similar  trend  to  the 
result  shown  in  Figure  7. 


Q- 


0  2  4  6  8  10  12  14  16  18  20 

SCR,  dB 

Figure  7:  Performance  of  amplitude  detection  and  log  detection  in  comparison  to  the  optimal 
detection,  i.e.,  the  intensity  detection. 


3.2  Clutter  with  Temporal  Correlation 

Definitions  of  temporal  and  spatial  correlations  as  well  as  their  calculations  for  the 
compound  K-distributed  clutter  are  discussed  in  Appendix  C.  We  begin  the  discussion  for 
temporal  correlation  in  this  subsection  and  leave  the  discussion  for  the  spatial  correlation 
to  Section  4.  When  measure  x  is  temporally  correlated,  the  associated  marginal  pdf 
becomes. 


(12) 


(13) 


where  M  f  is  the  normalised  covariance  matrix  of  clutter  with  respect  to  pulse,  the  matrix 
can  also  be  viewed  as  the  covariance  matrix  of  the  fast-varying  component  of  the  clutter 
(See  Appendix  C  for  details),  and  Ms  is  the  normalised  covariance  matrix  of  target  signal 
with  respect  to  pulse. 
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Following  the  derivation  process  shown  in  the  previous  subsection,  the  corresponding 
LTR  for  the  model  presented  by  (12)  and  (13)  is, 

2(x)  (14) 

It  can  be  denoted  by, 

/l(x)  =  —  xHs  (15) 

N 

where 

s  =  M_1x  (16) 

M-1  -(MMf+a2MsY  (17) 

The  diagram  of  such  an  optimal  detector  is  shown  in  Figure  8.  The  N-P  detector  correlates 
the  received  signal  x[n\  with  an  estimate  of  the  signal  s\ri\ .  It  is  therefore  termed  an 

estimate-correlator  (Kay,  1998,  Chapter  5),  and  s  is  usually  referred  to  as  a  Wiener  filter 
estimator  of  the  signal. 


Figure  8:  Estimate-correlator  for  non-coherent  detection  of  Gaussian  random  signal  in  compound 
K-distributed  clutter. 


3.2.1  Case  I:  Mf  =Ms 

If  both  clutter  and  target  signal  have  the  same  normalised  covariance  matrix, 
M ,  =  Ms  =  M0,  the  above  Wiener  filter  can  be  simplified  to, 

M_1  =  Mq1  (18) 
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This  is  the  well-known  whitening  processing.  The  coloured  (correlated)  measures  become 
white  (uncorrelated)  after  the  whitening  processing.  Since  the  data  becomes  uncorrelated 
after  the  whitening  processing,  the  detection  performance,  including  the  threshold,  false- 
alarm  rate  and  probability  of  detection  are  all  the  same  as  those  discussed  in  the  previous 
subsection. 


3.2.2  Case  II:  Mf  *  M, 


Unlike  clutter  data,  whose  parameters  can  be  measured  or  estimated  from  the  secondary 
data,  target  signals  are  normally  unknown.  If  there  is  no  prior  knowledge,  target  signals 
may  be  assumed  to  be  uncorrelated,  giving. 


M"1  =  {ju  M;  )-1  -  (ju  My  +  cr2l)-1 
The  matrix  inversion  lemma  leads  to. 


ivr1 


=  —  M31 
A 


f  i  i 

-  M/+  —  I 

VA  o- 


Y 


J 


—  M) 

A 


(19) 


(20) 


It  can  be  seen  from  (20)  that  owing  to  different  correlation  properties  of  clutter  and  target 
signals,  the  optimal  processor  cannot  form  a  filter  that  de-correlates  correlated 
measurements  into  fully  uncorrelated  measurements  under  either  H0  or  //  hypothesis. 
Instead,  the  processor  compromises  the  different  correlation  properties  of  clutter  and 
target  signal,  and  de-correlates  measurements  into  overall  least  correlation  to  achieve  the 
optimal  detection.  It  can  be  shown  that  M_1  ~  //"'Mi1  for  a  large  SCR  (strong  target)  case. 


Implementation  of  the  Wiener  filter  (19)  for  the  optimal  processor,  however,  encounters  a 
difficulty  as  it  requires  knowledge  of  the  target  signal's  intensity.  Because  there  is  no  prior 
knowledge  of  target  signals,  such  an  optimal  detector  is  difficult  to  implement. 

The  generalised  likelihood  ratio  test  (GLRT)  allows  a~  to  be  estimated  by  the  maximised 
likelihood  estimate  (MLE)  method.  The  MLE  of  a 2  is  obtained  by  maximising  the 
marginal  pdf  of  max  {in  p(x;  r,  //, )} . 


In  p(x;  t,H1)  =  -N  1n(2;r)  -  In  det(Mc,  +  <r2I  )-x"(m  c+<72i)-'x 


N  bi^-Iln^+rr2)-!^ 

n= 1 


Efxf 


^  n  2 

n= 1  /t  +  (7 


(21) 
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where  Mc  =  rMf/  An  and  Ej; ,  n  =  !,•••,  N ,  are  the  eigenvalues  and  eigenvectors  of  Mc. 
To  find  the  MLE  of  a2 ,  we  must  minimise  J(cj2 )  with  respect  to  <x2,  where 


Act2) 


N 

I 

n=\ 


ln(A„  +cr2)- 


EfV 


A, 


a 


(22) 


Differentiating  J((J2)  leads  to  a  non-linear  equation  of  a2,  and  there  is  no  general 
solution.  Therefore,  the  GLRT  detector  cannot  be  found  analytically. 

We  therefore  need  to  find  suboptimal  detectors.  One  of  suboptimal  detector  is  to  use 
a  M-/  (a  is  an  arbitrary  real  number),  because  the  filter  fully  de-correlates  the  clutter  to 
achieve  the  lowest  threshold  for  a  given  false-alarm  rate.  This  suboptimal  detector  also 
becomes  optimal  for  the  strong  target  case.  Selecting  a  =  1 ,  the  data  after  the  whitening 
processing  will  become  statistically  identical  to  the  case  of  uncorrelated  clutter  discussed 
in  Subsection  3.1  under  hypothesis  H0 ,  so  that  the  threshold  setting  will  be  unchanged. 
The  whitening  filter,  however  will  affect  the  measurement  of  target  signal  which  in  turn 
alters  the  detection. 


For  uncorrelated  Gaussian  target  signal  vector  xt,  the  mean  measurement  of  multi-looks  is 
E  {xf  x, }  =  N a 2 .  When  the  whitening  filter  is  used,  the  measurement  becomes. 


where  tr{)  denotes  trace  of  matrix  and  Xn ,  n  =  !,■■■  ,N ,  are  eigenvalues  of  IY1  ,  .  Because 
fr(lVI , )  =  X,)  i  \  =  N ,  we  immediately  have. 


N  1 

T—>N  (24) 

n=1  K 


The  equal  sign  holds  if  and  only  if  ^  =  •  •  •  =  XN  =  1  that  is  for  the  uncorrelated  case  of 
M  f  =  I .  Therefore,  when  clutter  is  correlated  and  the  whitening  filter  is  used  to  fully  de¬ 
correlate  the  correlated  clutter,  the  whitening  filter  over  the  uncorrelated  Gaussian  target 
signal  always  results  in  a  target  gain.  For  instance,  if  the  correlation  coefficient  is  a 
geometric  series,  i.e.,  pn=  p"  ,  n  =  0,---,N  —  l  (such  as  pn  =  exp(-/f  n)  ,  J3  >  0  is  a 
constant)  then  we  can  show, 

*to)-N+Zi??pf  {25) 
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The  corresponding  target  signal  gain  will  be. 


N  +  (N-2)\p 
A(l-|p|2) 


(26) 


The  proof  of  (25)  may  be  found  in  Appendix  D.  Therefore,  the  higher  the  correlation,  the 
higher  the  target  gain.  The  pdf  of  z  =  xf  M~'x(  will  not  be  gamma,  and  can  be  derived 
using  the  method  given  in  Appendix  E  if  desired. 

The  reason  we  obtain  a  target  signal  gain  by  applying  the  whitening  filter  may  be 
explained  in  this  way.  When  both  target  signal  and  clutter  are  uncorrelated,  the  intensity 
estimation  is  the  only  way  to  detect  the  presence  of  target  signal.  However  when  clutter  is 
correlated,  its  spectrum  is  not  uniformly  distributed  anymore,  whereas  the  uncorrelated 
target  signal  has  a  uniformly  distributed  spectrum.  In  other  words,  the  difference  in 
spectra,  if  being  utilised  can  improve  the  detection  (or  equivalently,  providing  a  target 
signal  gain).  The  function  of  the  whitening  filter  in  the  spectral  domain  is  to  multiply  a 
least  coefficient  to  the  largest  power  spectral  component  and  a  largest  coefficient  to  the 
least  power  spectral  component  so  that  the  filtered  data  has  a  uniformly  distributed  power 
spectrum.  Applying  the  same  coefficients  to  the  uniformly  distributed  power  spectral 
components  (uncorrelated  target  signal),  however,  provides  a  gain.  As  an  illustration. 
Figure  9  depicts  power  spectra  of  the  correlated  clutter  and  the  uncorrelated  target  signal, 
respectively  ( N  =  6).  As  shown  in  the  figure,  the  original  SCR  is. 


SCR 


l+l+l+l+l+l 

2  +  5  +  3  +  1  +  0.5  +  0.5 


(27) 


The  whitening  processing  is  to  find  a  set  of  coefficients  so  that  the  spectrum  of  clutter  after 
filtering  becomes  uniformly  distributed  while  maintaining  the  integral  of  the  power 
spectrum  (the  total  power)  unchanged.  Such  a  set  of  coefficients  is, 
{l  2/5  2/3  2  4  4}.  However  the  same  set  coefficients  when  applied  to  the 
spectrum  of  uncorrelated  target  signal  generates  a  gain,  as  the  SCR  after  filtering 
processing  becomes. 


SCR  = 


1  +  2/5  +  2/3  +  2  +  4  +  4 


(28) 


Comparing  (27)  and  (28),  the  generated  gain  is  2  (3dB).  From  this  example,  we  can  also  see 
that  the  higher  the  correlation,  the  higher  the  target  signal  gain,  as  indicated  previously. 
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(b) 


Figure  9:  Power  spectra  of  (a)  correlated  clutter  and  (b)  uncorrelated  target  signal. 


In  conclusion,  if  target  signal  and  clutter  have  different  correlation  properties,  and  there  is 
no  prior  knowledge  about  the  target  signal,  the  GLRT  detector  cannot  be  found.  One  of  the 
sub-optimal  detectors  for  the  uncorrelated  Gaussian  target  signal  embedded  in  the 
correlated  K-distributed  clutter  is  to  use  the  inverse  of  the  normalised  covariance  matrix  of 
clutter  as  the  whitening  filter  (this  sub-optimal  detector  becomes  optimal  for  large  SCR 
targets).  After  the  whitening  processing,  the  associated  threshold  for  a  given  false-alarm 
rate  will  be  the  same  as  that  for  the  uncorrelated  clutter.  On  the  other  hand,  the  whitening 
process  provides  a  target  signal  gain  (the  value  of  the  gain  is  dependent  on  the  correlation 
properties  of  the  clutter),  and  hence  results  in  a  higher  probability  of  detection  compared 
to  the  same  target  signal  embedded  in  uncorrelated  clutter.  In  another  words,  the  ROC 
curves  discussed  in  Subsection  3.1  is  for  the  worst  case  scenario. 

To  demonstrate.  Figure  10  compares  the  performance  of  the  sub-optimal  detector  using 
m:1  as  the  whitening  filter  for  detecting  uncorrelated  Gaussian  target  signals  embedded 

in  correlated  K-distributed  clutter.  The  blue  ROC  curve  is  analytically  calculated  for  the 
uncorrelated  Gaussian  target  signal  embedded  in  the  uncorrelated  K-distributed  clutter 
with  a  shape  parameter  of  v  =  1 .2 .  The  green  asterisks  are  the  result  of  the  Monte  Carlo 
simulation,  which  matches  the  theoretical  results.  The  second  ROC  curve  (red  line  with 
small  circles)  is  the  result  of  the  Monte  Carlo  simulation  for  a  case  of  the  uncorrelated 
Gaussian  target  signal  embedded  in  the  correlated  K-distributed  clutter  with  the 
correlation  coefficients  of  pn  =  exp(-/?),  n  =(),■•■,  N  —  I  .  The  third  ROC  curve  (broken 
purple  line  with  small  circles)  is  the  result  of  the  Monte  Carlo  simulation  for  a  case  of  the 
uncorrelated  Gaussian  target  signals  embedded  in  the  correlated  K-distributed  clutter  with 
the  correlation  coefficients  of  pn  =  exp(-«  /  2) ,  n  =(),■■•,  N  —  I .  As  illustrated  earlier,  since 
the  whitening  processing  produces  a  target  signal  gain  for  the  uncorrelated  target  signals, 
the  detection  for  the  two  correlated  clutter  cases  are  improved.  According  to  (26),  the 
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associated  target  signal  gains  are  1.28  (l.ldB)  and  2.05  (3.1dB)  for  the  two  correlated  cases, 
respectively.  Knowing  the  target  signal  gain,  the  associated  probability  of  detection  can 
also  be  approximately  found  from  the  probability  of  detection  for  the  uncorrelated  case2. 
For  instance,  the  target  gain  for  the  second  correlated  case  is  3.1dB,  which  means  the 
probability  of  detection  for  a  target  having  a  SCR  of  lOdB  for  the  correlated  case  will  be 
close  to  the  probability  of  detection  for  a  target  having  an  SCR  of  13.1dB  for  the 
uncorrelated  case,  i.e.,  the  third  ROC  curve  can  be  approximately  obtained  by  horizontally 
left-shifting  the  first  ROC  curve  by  3.1dB.  Similarly,  the  second  ROC  curve  can  be 
approximately  obtained  by  horizontally  left-shifting  the  first  ROC  curve  by  l.ldB. 


1 
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Figure  10:  ROC  comparison  between  uncorrelated  and  correlated  K-distributed  clutter  embedded 
with  uncorrelated  Gaussian  target  signals.  Two  correlated  cases  are  shoivn:  Case  I  with 
correlation  coefficients  of  pn  =  exp(-«) ,  n  =  0,  •  •  ■ ,  N  - 1,  and  Case  II  with  correlation 

coefficients  of  pn  =  exp(-«)  ,  n  =  0,---,N  -1  (Pfa=  10“6  ,  p  =  1  ,  v  =  1.2  and 
N  =  10). 


If  target  signal  is  temporally  correlated  and  has  a  normalised  correlation  covariance  matrix 
[VI .  that  is  known  a  priori,  the  target  signal  mean  using  M_l  as  the  whitening  filtering  can 

be  calculated  by, 

^M^x^M^x,^  R-(m 's/2M ’/2  )<t2  (29) 


2  The  pdf  of  the  target  signal  after  whitening  processing  slightly  differs  from  the  multi-look  gamma 
distribution,  resulting  in  a  slightly  different  detection. 
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The  corresponding  target  signal  gain  then  can  be  calculated. 


4.  Detection  Against  Spatially  Correlated  Clutter 

Spatial  correlation  is  governed  by  the  slowly-varying  component  of  the  compound 
K-distributed  clutter.  The  spatial  correlation  is  the  correlation  with  respect  to  range  bin, 
while  the  temporal  correlation  discussed  earlier  is  the  correlation  with  respect  to  pulse. 
Details  of  how  to  discriminate  the  spatial  correlation  from  the  temporal  correlation,  as  well 
as  how  to  calculate  /  estimate  the  correlation  coefficients  are  given  in  Appendix  C.  For 
simplicity,  the  discussion  in  this  section  only  considers  the  spatial  correlation,  as  the 
temporal  correlation  has  been  discussed. 

The  difference  between  Gaussian  clutter  and  the  compound  K-distributed  clutter  lies  that 
the  texture  (the  mean  of  the  slowly-varying  component)  is  a  constant  for  the  former  and 
fluctuates  for  the  latter.  If  the  local  texture  r  were  known,  the  corresponding  N-P  optimal 
detector  would  be  (when  there  is  no  temporal  correlation). 


^(x) 


1  1  N 

T  N  «=1 


<  r 

Ho 


(30) 


where  the  threshold  y  is  determined  by  the  false-alarm  rate,  and. 


1  fa 


1 

i\N) 


F  {N,Ny) 


(31) 


r(u,h)  is  the  incomplete  gamma  function3,  defined  as,  l’(u, b)  =  j  \ta  'e  ' dt .  The  optimal 

b 

detector  of  (30)  can  be  directly  derived  from  the  pdf  of  multi-look  uncorrelated  Gaussian 
clutter  given  by  (B5).  The  detector  (30)  is  referred  to  as  the  'ideal  CFAR'  (Watts  1985;  Ward 
et  al.  2006),  and  its  performance  is  identical  to  the  optimal  detector  discussed  in  Section  3 
for  the  Gaussian  clutter  case. 

In  the  previous  Section,  the  slowly-varying  component  (texture)  is  assumed  to  be  fully- 
correlated  (remain  constant)  during  the  multi-pulse  collection  for  the  same  range  bin,  but 
fully-uncorrelated  from  range  bin  to  range  bin.  When  the  local  texture  r  is  spatially 
uncorrelated,  its  best  estimate  is  its  global  mean  //  .  The  optimal  and  near-optimal 


3  Matlab  defines  the  incomplete  gamma  function  in  a  different  way,  as,  Y{b,  a) 
As  a  result,  Pfa  =  1  —  I  ’(/V  y,  N) . 


1  b 

- f  t^e-'dt. 

V(  a  )  J0 
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processors  presented  in  Section  3  implicitly  use  the  global  mean  as  the  estimate  of  the  local 
texture  r  because  it  is  the  best  estimate  when  there  is  no  correlation. 


4.1  Estimation  of  Local  Texture 

The  texture  of  sea  clutter  is  believed  to  be  highly  correlated  to  wave  and  swell  structures  of 
the  sea  (Ward  et  al.  2006,  Chapter  2),  and  some  waves  can  be  metres,  tens  or  hundreds  of 
metres  long.  This  leads  to  the  slowly-varying  component  to  be  possibly  spatially 
correlated.  For  spatially  correlated  sea  clutter,  the  local  mean  f  estimated  by  the 
neighbouring  range  bins  of  CUT  is  often  a  better  estimate  than  the  global  mean  fj. .  As  a 
consequence,  using  the  detector  (30)  with  the  estimate  f  as  the  replacement  of  the 
unknown  parameter  r  may  result  a  better  detection.  This  is  known  as  the  CFAR  gain 
(Watts  1996;  Watts  et  al.  2007;  Watts  1987;  Ward  et  al.  2006).  Whether  a  CFAR  gain  is 
achievable  depends  on  three  conditions: 

1.  Clutter  has  fluctuating  texture  (i.e.,  compound  non-Gaussian  distributed); 

2.  The  texture  is  spatially  correlated;  and 

3.  The  local  texture  is  estimated  correctly. 

The  CFAR  gain  in  turn  depends  on  the  accuracy  of  the  estimate  of  the  local  mean  f ,  and 
the  spikiness  of  the  clutter.  If  clutter  is  spiky,  and  its  spatial  correlation  is  high,  obtaining  a 
CFAR  gain  is  possible  by  a  proper  estimation  for  the  local  texture. 

It  should  be  pointed  out  that  the  application  of  the  above  method  is  not  limited  to  the  case 
of  compound  K-distribution.  In  fact,  any  compound  non-Gaussian  distribution  whose 
texture  may  have  different  pdfs  but  possess  the  similar  characteristic  (e.g.,  constant  during 
the  multi-pulse  collection  and  fluctuating  from  range  bin  to  range  bin  and  spatially 
correlated)  falls  in  this  category. 

As  stated  above,  the  better  the  estimate  of  the  local  texture  f ,  the  better  the  performance 
of  the  detector.  Therefore,  finding  the  ways  that  make  better  estimates  of  the  local  texture 
t  improves  the  performance  of  the  detector. 

Following  the  idea  given  by  Bucciarelli  et  al  (1996),  we  use  the  linear  autoregressive  (AR) 
technique  to  estimate  the  local  texture  f  .  Because  E{z}  =  E{r}E{\  xf  |2}  =  E{t}  ,  we  can 
use  the  clutter  intensity  z  to  estimate  the  underlying  texture  r  . 
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(a)  Before  multi-look  averaging  processing  (b)  After  multi-look  processing 

Figure  11:  Use  neighbouring  measures  to  estimate  CUT  (the  case  of  multi-look  in  azimuth). 


Figure  11  depicts  the  use  of  neighbouring  cells  to  estimate  the  texture.  As  shown  in  Figure 
11  (b),  for  a  multi-look  CUT  z[0]  =  J  x0[/]  |2,  we  want  to  use  its  neighbouring  cells 

to  estimate  its  value.  Using  the  AR  technique,  the  estimate  of  z[0]  may  be  expressed  as  a 
linear  combination  of  its  neighbouring  cells,  as  (for  a  general  stationary  and  symmetrical 
random  process), 

z[0]  =  X  wn  (z[n]  +  z[-n ])  (32) 

n=l 

where  L  is  the  length  of  the  one-side  estimation  window.  The  unknown  parameters,  wn , 
n  =  !,•••, L ,  should  satisfy, 

min  (|z[0]-z[0]|  \  (33) 

wn,n=l,---,L  \  / 

Among  the  L  unknown  parameters,  wn  ,  n  =  1,  •  •  • ,  L  ,  however,  only  L  —  1  are 
independent,  as  (z[0]  )  =  £  w„  ((z[«]  +  z[-«]))  leads  to, 

n= 1 

(34) 

n= 1  Z 

Without  loss  of  generality,  we  rewrite  (34)  as, 

wi  =  A  “  2  w„  (35) 

2  n= 2 

The  minimisation  of  (33)  can  be  found  by  use  of  the  Lagrange  theorem,  as. 
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Inserting  (32)  and  (34)  into  (36),  and  after  some  manipulation,  we  have. 


\(2rlk~  2rlx  ~  rjk-i  ~  7lk+\  +  *12  +  r/o ) 

=  Z  wn  b„+k  +  %-k\  -  r/n-x  -  nn+x  -  Vk-x  -  nk+l  +r/2+n0) 


k  =  2,--,L 


(36) 


(37) 


where  r)k  = 


E{z[i]z[i  +  k]}-E2{z} 
var(z) 


,  k  =  0, 1,  •  •  ■ ,  is  the  correlation  coefficient  of 


z  . 


Changing  the  value  of  k  gives  L—  1  linear  equations,  which  provides  a  unique  solution 
for  wn ,  n  =  2,---,L  .  Together  with  (35),  the  all  unknown  parameters,  wn ,  n  =  1,2,---,  A, 
can  be  uniquely  obtained. 

It  can  be  shown  that  for  spatially  uncorrelated  data,  the  weights  become  equal  and 
wn  =1/(2 L) ,  n=l,---,L  .  That  is  consistent  with  the  well-known  CA  processing.  If, 
however,  clutter  is  spatially  correlated,  the  optimal  weights  are  unequal. 


Considering  extended  target  signals  which  may  occupy  a  few  range  cells,  radar  engineers 
often  prefer  to  exclude  a  few  guardian  cells  next  to  CUT  when  estimating  CUT  using 
neighbouring  cells,  which  is  shown  in  Figure  12.  The  corresponding  optimal  weights  can 
be  found  accordingly  using  the  following  linear  equations. 


CUT 


d-L-g] 

d-i-gl 

z[0] 

z[l+g] 

z[L+g] 

Range 


\  X 

Guardians 


Range 


Figure  12:  CUT  estimation  using  neighbouring  cells  with  guardian  cells. 


\  (2?7r+g  -  2flx+g  -  Vk- 1  -  +  th+2g  +  %) 

=  X  K  b„+k+2g  +  %-k\  -  Vn-X  -  Vn+  l+2g  -  Vk-X  ~  V k+X+2g  +  V2+2g  +  V 0  ) 

n= 2 


k  =  2,---,L  (38) 


where  g  is  the  number  of  guardian  cells  at  each  side  excluded  in  the  estimation 
processing. 
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4.2  Results 

After  the  local  texture  of  CUT  is  estimated,  we  can  now  use  the  suboptimal  detector. 


T(x) 


I  1  »  , 

--H\x[n]  |2 
r0  N  n= i 


H\ 

<  Y 
H0 


(39) 


to  carry  out  the  detection,  where  f0  =  z[0]  is  the  estimate  of  CUT's  texture  obtained  by  the 
previously  subsection.  However,  since  the  pdf  of  T(x)  is  unknown,  the  threshold  for  the 
above  suboptimal  detector  is  also  unknown,  and  has  to  be  simulated  using  the  Monte 
Carlo  simulation. 

Using  the  Monte  Carlo  simulation,  the  performance  of  the  detector  (39)  is  determined  by 
the  following  two  steps. 

1.  Generating  spatially  correlated  K-distributed  clutter  data  accordingly,  estimating 
its  local  texture,  calculating  the  test  ratio  using  (39),  and  numerically  determining 
the  threshold.  The  simulation  used  AxlO9  samples  to  determine  the  threshold, 
which  ensures  an  absolute  relative  error  less  than  5%  for  95%  of  time  for  a  false- 
alar  m  rate  of  10“'’  ( Pfa  =  10“6 )  (Kay  1998,  Chapter  2). 

2.  Generating  clutter  data,  adding  target  signal  data  and  determining  the  probability 
of  detection  accordingly.  The  detection  calculation  used  N  x  106  samples. 

The  detection  improvement  (CFAR  gains  achieved)  varies  depending  on  the  correlation. 
Below  we  show  two  numerical  samples  to  demonstrate  the  improvement  of  using  the 
optimal  weights  to  estimate  the  local  texture.  We  assume  there  is  no  temporal  correlation 
because  such  correlation  can  be  decorrelated  using  the  technique  described  in  Section  3. 
The  spatial  correlation  of  the  K-distributed  clutter  is  resulted  from  the  correlation  of 
gamma  distributed  texture. 

In  the  first  sample,  we  assume  the  gamma  distributed  clutter  texture  has  a  correlation 
coefficient  of, 

Pk  =[0.7  +  0.3cos(0.12^A:)]e“i/!2  it  =  0,1,-  (40) 

where  k  denotes  the  number  of  lagged  range  bins.  After  the  modulation  of  the  texture 
with  the  fast-varying  Gaussian  component,  and  multi-look  (independent  looks  in  azimuth) 
intensity  averaging  processing,  the  resultant  multi-look  K-distributed  clutter  has  a 
correlation  of  (refer  to  Appendix  C  for  details). 
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=  1  and  77,  =  pk - — —  k=  1,2—  (41) 

v  +  N +  1 

Values  of  v  =  1.2  and  N  =  10  was  used  in  the  Monte  Carlo  simulation.  The  texture 
estimated  by  the  optimal  weights  determined  by  (37)  is  much  better  than  the  estimate  of 
the  CA  window  when  clutter  is  spatially  correlated.  Figure  13  shows  the  estimates  of  CUT 
using  these  two  methods  in  comparison  with  the  true  value.  It  can  be  seen  that  the 
estimates  by  the  optimal  weights  follow  the  fluctuation  much  better  whereas  the  estimates 
of  the  CA  processing  do  not  seem  to  follow  the  fluctuation  closely. 


Figure  13:  CUT  estimate  by  use  of  the  optimal  weights  and  the  CA  window,  respectively.  Data 
parameters  are :  v  =  1 . 2 ,  jj.-  1.0,  N  =  10,  and  L  =  8 .  The  spatial  correlation  is  given 
by  (40)  and  (41). 


Using  the  optimal  estimates  of  the  texture  for  the  suboptimal  detector  (39),  however,  did 
not  show  a  significant  improvement  of  the  detection.  A  close  examination  indicated  that 
small  values  of  the  estimate  f  could  lead  to  large  values  of  A(x)  (see  (39)),  which,  in  turn, 
forces  the  threshold  to  be  set  high  to  maintain  the  desired  false-alarm  rate.  We  recall  that 
false-alarms  are  often  caused  by  sea  spikes,  e.g.,  the  high  clutter  returns,  if  the  detector 
uses  a  fixed  threshold,  and  the  low  returns  are  not  the  concern.  Therefore,  it  seems  that 
there  is  no  need  for  a  close  estimate  for  low  clutter  returns,  and  an  accurate  estimate  is 
required  only  for  high  clutter  returns. 

To  improve  the  detector's  performance  The  estimate  f  is  modified  as, 

i  -p[°] 

0  \ju/2  if  z[0]  <  /u!  2  1  ’ 
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The  profile  of  the  modified  estimates  of  the  clutter  texture  is  shown  in  Figure  14.  The 
modification  effectively  eases  the  false-alarms  produced  by  low  clutter  returns,  and  hence 
greatly  improves  the  detection.  Our  simulation  indicated  that  the  overall  performance  is 
not  very  sensitive  to  the  selection  of  the  cut-off  of  f  (for  instance,  (42)  selects  jul  2  as  the 
cut-off  value).  In  general,  a  lower  cut-off  value  improves  detection  a  little  for  low  SCR 
targets  at  a  small  sacrifice  of  the  detection  of  high  SCR  targets  whereas  a  higher  cut-off 
value  improves  detection  a  little  for  mid  and  high  SCR  targets  at  a  sacrifice  of  the  detection 
of  low  SCR  targets. 


Figure  14:  Profile  of  the  modified  texture. 


The  performance  of  the  suboptimal  detector  (39)  using  the  modified  texture  estimates 
given  by  (42)  is  shown  in  Figure  15.  Compared  to  the  fixed  threshold  (which  assumes  the 
mean  of  clutter  be  known),  the  suboptimal  detector  greatly  improves  the  detection  (the 
improvement  is  referred  to  as  the  CFAR  gain  by  Watts  and  others  (Watts  1996;  Watts  et  al. 
2007;  Watts  1987))  and  provides  a  CFAR  gain  about  6dB  at  Pd  =  0.5 .  The  detection  using 
CA  processing  also  achieves  some  CFAR  gain  (about  2dB  at  Pd  =  0.5),  but  the  proposed 
detector  performs  far  better. 
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Figure  15:  ROC  comparison  between  the  fixed-threshold  and  the  adaptive  thresholds  determined  by 
the  optimal  zveights  and  a  CA  window  with  L  =  8  for  the  spatially  correlated  K- 
distributed  clutter  with  the  correlation  given  by  (40)  and  (41). 


In  the  second  numerical  example,  we  assume  the  gamma  distributed  clutter  texture  has  a 
correlation  coefficient  of, 

pk=e~kn  A:  =  0,1, —  (43) 

Obviously,  this  correlation  is  much  lower  than  that  of  the  first  example.  Again,  the  multi¬ 
look  K-distributed  clutter  will  have  a  correlation  coefficient  specified  by  (41).  The 
performance  of  the  suboptimal  detector  (39)  using  the  modified  texture  estimate  given  by 
(42)  is  shown  in  Figure  16,  together  with  the  performances  of  the  CA  window  and  the 
fixed  threshold.  It  can  be  seen  that  the  proposed  processor  is  still  able  to  provide  a 
moderate  CFAR  gain  (about  2dB  at  the  Pd  =0.5),  whereas  the  CA  processor  performs 
approximately  the  same  as  the  fixed  threshold. 
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Signal-to-clutter  ratio  (dB) 


Figure  16:  ROC  comparison  between  the  fixed-threshold  and  the  adaptive  thresholds  determined  by 
the  optimal  weights  and  a  CA  window  with  L  =  8  for  the  spatially  correlated  Re¬ 
distributed  clutter  with  the  correlation  given  by  (43)  and  (41). 


It  is  worth  noting  that  the  achievable  CFAR  gain  depends  on  the  correlation  and  the  shape 
parameter.  The  higher  the  correlation  and  smaller  shape  parameter,  the  greater  the  CFAR 
gain.  If  there  is  no  correlation,  there  will  be  no  CFAR  gain  at  all.  Instead,  a  CFAR  loss  will 
occur.  Because  for  spatially  uncorrelated  clutter,  the  best  estimate  of  the  local  mean  is  the 
global  mean,  and  using  a  limited  size  of  the  CA  window  to  estimate  the  local  mean 
inherently  associates  with  a  loss  (the  optimal  weights  becomes  identical  to  the  CA  window 
for  the  spatially  uncorrelated  clutter). 

'Ideal  CFAR'  proposed  by  Watts  (1985)  is  the  theoretical  upper  bound  a  CFAR  processor 
could  achieve  for  spatially  correlated  non-Gaussian  clutter.  However,  how  to  find  the 
unknown  'true'  clutter  texture  remains  unsolved.  This  section  uses  the  linear  AR  technique 
to  estimate  the  unknown  clutter  texture.  Theoretically,  the  estimate  is  the  optimal  for  the 
given  size  of  the  window  and  the  linear  regression  model.  The  estimated  texture  is  then 
modified  to  mitigate  possible  false-alarms  caused  by  low  returns.  It  has  shown  through 
numerical  examples,  the  proposed  suboptimal  detector  is  able  to  provide  a  much  higher 
CFAR  gain  compared  to  the  traditional  CA  processor  for  spatially  correlated  non-Gaussian 
clutter. 
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5.  Summary 

Optimal  non-coherent  radar  detection  of  Gaussian  targets  (Swerling  II  model)  embedded 
in  the  compound  K-distributed  clutter  has  been  investigated.  The  derived  optimal 
detector,  under  the  sense  of  the  Neyman-Pearson  principle  is  the  well-known  square-law 
detector.  This  is  because  the  underlying  texture  of  the  compound  K-clutter  remains 
unchanged  during  the  multi-pulse  averaging  processing  period  according  to  the 
assumption,  so  that  the  clutter  actually  undergoes  a  Gaussian  random  process  during  the 
period  of  multi-pulse  collection.  While  the  derived  optimal  detector  is  not  new,  the 
derivation  itself  has  a  guiding  meaning.  As  the  detector  has  been  rigorously  derived  for 
the  first  time  using  the  Neyman-Pearson  principle,  it  means  there  does  not  exist  any 
detector  that  would  perform  better  for  the  given  conditions.  Other  detectors,  such  as 
multi-look  amplitude  detector,  multi-look  log  detector  and  the  like  are  inherently 
associated  with  some  detection  loss  compared  to  the  optimal  detector,  the  multi-look 
intensity  detector. 

A  sub-optimal  detector  (it  becomes  optimal  for  large  signal-to-clutter  ratio  targets)  has 
been  proposed  for  temporally  correlated  K-distributed  clutter.  The  normalised  covariance 
matrix  of  clutter  is  used  to  de-correlate  the  clutter.  The  whitening  processing  provides  a 
target  signal  gain,  resulting  in  an  improved  detection.  The  higher  the  correlation,  the 
larger  the  target  signal  gain.  The  target  gain  results  from  the  difference  in  spectra  of 
correlated  clutter  and  the  uncorrelated  target  signals  (a  non-uniform  spectrum  of  clutter 
against  a  uniform  spectrum  of  target)  that  provides  a  second  characteristic  (in  addition  to 
the  intensity)  for  discriminating  the  target  from  the  clutter.  The  occurrence  of  uncorrelated 
Gaussian  targets  embedded  in  uncorrelated  compound  clutter  represents  the  worst 
scenario  in  terms  of  detection. 

If  the  clutter  is  spatially  correlated,  using  a  limited  number  of  neighbouring  range  bins  to 
estimate  the  local  texture  can  provide  a  CFAR  gain  and  improve  the  detection.  This  report 
has  proposed  the  use  of  the  linear  AR  technique  and  derived  the  optimal  weights  for 
estimating  the  local  texture  of  clutter.  The  AR  estimation  is  optimal  under  the  linear 
assumption  and  results  in  better  results  than  the  cell-averaging  estimation.  The  AR 
estimates  are  further  modified  to  mitigate  possible  false-alarms  caused  by  low  clutter 
returns.  It  has  shown  that  the  proposed  processor  greatly  improves  the  detection 
compared  to  the  traditional  CA  processor  for  spatially  correlated  clutter.  The  higher  the 
correlation,  the  higher  the  CFAR  gain.  However,  when  clutter  is  spatially  uncorrelated,  the 
global  mean  is  the  best  estimate  of  the  local  texture. 

A  summation  of  the  optimal  or  near  optimal  non-coherent  process  for  detecting  Gaussian 
targets  embedded  in  K-distributed  clutter  that  may  be  temporally  and  spatially  correlated 
is  provided  in  Appendix  A  as  a  quick  reference  for  radar  engineers  implementing  the 
detection  process. 
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Appendix  A:  Summation  of  Optimal  or  near  Optimal 
Non-coherent  Detection  Process 


This  appendix  is  a  summation  of  the  optimal  or  near  optimal  non-coherent  radar  process 
discussed  in  this  report.  It  serves  as  a  quick  reference  for  radar  engineers  implementing 
the  detection  process. 

Assumptions: 

1.  Clutter  has  a  compound  K-distribution  (thermal  noise  is  included),  and  may  have 
temporal  correlation  (with  respect  to  pulse)  and  spatial  correlation  (with  respect  to 
range).  Parameters  of  the  K-distribution  as  well  as  the  correlations  are  known  or 
have  been  estimated  using  the  sample  data. 

2.  Target  is  point-like  and  has  a  Gaussian  distribution. 

For  each  cell  under  test  (CUT),  its  neighbouring  cells  under  the  consideration  is  depicted 
in  Figure  Al.  The  neighbouring  cells  under  the  consideration  are  assumed  to  be  target- 
free. 
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Figure  Al:  Using  neighbouring  cells  to  estimate  CUT  (the  case  of  multi-look  in  pulse  and  multi¬ 
look  in  range). 


First,  the  temporal  correlation  is  decorrelated  and  cell  values  are  averaged  for  each  range 
bin.  This  process  is  called  multi-look  process. 

z;=^xf  Mx,  /  =  ±l,-,±£  (Al) 

where 
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x,=[x,[l]  •••  X/[jV]f  (A2) 

The  calculation  of  M_l  is  given  by  (20).  If  target  signal  is  uncorrelated  or  unknown,  M_l 
becomes, 

M”1  =  /T'lVI}1  (A3) 

where  fi  is  the  clutter  mean  and  IN'I ,  is  the  normalised  temporal  correlation  matrix  of 
clutter.  If  clutter  is  temporally  uncorrelated,  (Al)  becomes, 

zi=\ -XhM’  l  =  ±l,-",±L  (A4) 

n= 1 

The  above  process  reduces  the  two-dimensional  data  window  to  a  one-dimensional  data 
window  as  shown  in  Figure  A2. 
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Figure  A2:  Data  window  after  multi-look  processing. 


The  second  stage  is  to  estimate  z[0]  by  symmetrical  weights,  as, 
z[0]  =  YJW„{z[n]  +  z[-n]) 

n= 1 


The  weights  are  calculated  by. 
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Where 

K  =  j  {2tlk+g  -  2Vi+g  -  Vk-\  -  r/k+i+2g  +  n2+2g  +n0)  k  =  2,  ■  ■  ■ ,  L 


(A5) 


(A6) 


(A7) 
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a,„  = 


Vn+k+2g  +  %-k\  -  Vn-\  -  >7,i+l+2g  “  Vk-l  ~  H  k+l+2g  +  H  2+  2g  +V0  k,tl  =2  ,—,L 


(A8) 


,  E{z[i]z[i  +  k]\-E2{z)  ,  .  ,  ,  .  ....  , 

where  rjk  =  — - - - - — ,  k  =  (J,  1,  •  •  •  is  the  spatial  correlation  coethcient  ot  z,  and 

var(z) 

g  is  the  number  of  guardian  bins  on  each  of  side  of  CUT.  The  correlation  is  assumed  to  be 
symmetrical  with  respect  to  CUT. 

Finally  we  calculate  f0  by. 


Tn  =  ■ 


\z[0]  if  z[0]  >  /j,  1 2 
[ju/2  ifz[0]<///2 


(A9) 


The  detection  is  then  performed  simply  by  comparing  f0/ ju  to  a  threshold  y  that  is 
determined  by  the  false-alarm  rate  to  declare  absence  or  presence  of  target  in  the  CUT.  The 
whole  detection  process  is  depicted  in  Figure  A3. 


Figure  A3:  Block  diagram  of  optimal  or  near  optimal  non-coherent  detection  of  Gaussian 
targets  embedded  in  compound  K-distributed  clutter. 


If  we  let  wn  =  ,  n  =  !,•••, L ,  then  the  spatial  processing  becomes  cell-average  CFAR 


(CA-CFAR)  in  range.  If  there  is  no  spatial  correlation  (correlation  in  range),  there  is  no 
need  to  estimate  z[0]  using  multiple  range  cells.  The  best  estimation  of  z[0]  is  ju  in  this 
case. 
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Appendix  B:  Approximation  of  Gaussian  Signal  Added 

in  K-distributed  Clutter 


The  compound  K-distribution  is  assumed  to  have  two  components,  a  fast-varying 
component,  called  speckle,  that  is  a  zero  mean,  unit  variance  complex  Gaussian  variable, 
modulated  by  a  slowly-varying  component,  called  the  underlying  mean  or  texture  whose 
intensity  is  gamma  distributed.  The  intensity  of  such  K-distributed  clutter  has  a  pdf  of. 


P(z)  =  ]p(z\T)p(T)dT  (Bl) 

o 


where  z  =|  xc  |2  is  the  intensity  of  clutter,  p(z  |  r)  is  the  marginal  pdf  of  z  for  a  given 
underlying  mean  r  and  p(r)  is  the  pdf  of  the  underlying  mean.  Accordingly, 


p{z  |  t)  =  —  exp(-z  /  r) 
r 

(B2) 

P(T)  =  r,  exP(  bT ) 

r(v) 

(B3) 

v  2 

where  v  is  the  shape  parameter,  b= — ,  ju  =  E{  |  xc  |  }  is  the  mean  clutter  intensity. 

M 

Inserting  (B2)  and  (B3)  into  (Bl),  we  have  the  well-known  K-distribution  of 

p(z)  =  "FTt  (bzT'~lU2  Kv-1  (2Vfe )  (B4) 

nv> 

For  multi-look  processing,  if  during  the  multi-look  sampling  period,  the  slowly-varying 
component  remains  unchanged,  the  resultant  multi-look  processing  is  simplified  to  a 
multi-look  processing  of  fast-varying  component  modulated  by  the  unchanged  gamma 
distribution. 

If  the  fast-varying  component  varies  independently,  i.e.,  the  samples  are  uncorrelated,  the 
mean  of  N  exponentially  distributed  samples  is  a  gamma  distribution,  as. 

Inserting  (B5)  and  (B3)  into  (Bl),  one  obtains  the  N-look  K-distribution, 
p(z) =  nw)2r(>,)z  (Nbzftni K,_„ (2 4Wz )  (B6) 


■  /  AT 


(B5) 
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1  Jv 

where  z  =  —  X  z[h]  •  When  N  =  1 ,  (6)  simplifies  to  (B4). 
N  n=l 


Now  consider  the  distribution  of  a  combination  of  thermal  noise  (or  in  general  a  Gaussian 
signal)  and  K-distributed  clutter.  Assuming  the  Gaussian  signal,  and  clutter  each  is 
independently  random  process,  we  denote  Gaussian  signal  Xt~CN{ 0,  <r2),  and  the 
conditional  clutter  Xc  \  r  ~  CN (0,  r) .  Since  they  are  independent  processes,  and  Xt  is 
independent  of  r ,  we  have, 

(Xt+Xc)\r~CN(0,  t  +  ct2)  (B7) 


or. 


p(z 1 1  y) 


l 

y  +  <r2 


f 

exp 

v 


zi 

y  +  <t2 


\ 

) 


(B8) 


where  zx  =  xt  +  xc 


is  the  intensity.  The  pdf  of  is  given  by. 


oo 


P{Zl)  =  \p(zx  \T)p(T)dT 

0 


(B9) 


Unfortunately,  the  above  integral  does  not  have  a  closed-form  solution.  Therefore,  the 
distribution  of  the  Gaussian  target  signals  (or  Gaussian  noise)  embedded  in  the  K- 
distributed  clutter  is  unknown.  For  extreme  cases  we  note  that  the  distribution  will  still  be 
the  K-distribution  for  a  very  low  SCR  case  and  exponential  for  a  high  SCR  case.  The  latter 
is  also  a  K-distribution  with  the  shape  parameter  of  infinity.  Watts  (1987)  has  studied  the 
distribution  of  the  combined  K-distributed  clutter  plus  Gaussian  thermal  noise.  He  has 
approximated  the  combined  distribution  as  a  new  K-distribution,  with  a  modified  shape 
parameter  and  a  modified  mean.  Following  this  idea,  we  approximate  the  distribution  of 
(B9)  to  a  K-distribution  with  new  parameters  of  (jUt ,  v[ )  whose  values  are  derived  using 
moment  methods  in  the  following. 

It  is  known  that  for  the  K-distribution  given  by  (B4),  the  shape  parameter  is  a  function  of 
the  normalised  moments.  For  instance,  for  the  second  normalised  moment, 


2(\\  +D  _  E{z2} 
vi  E2{zx} 


(BIO) 


where  zx=\xt+xcf  ,  once  E\z]\  and  E{z2}  are  found,  the  corresponding  v,  can  be 
derived  accordingly,  as  the  distribution  of  zx  is  assumed  to  be  K. 

E{z i)  =  E{ |  x,  +  xr  |2}  =  E{ |  x,  |2}  +  E{ |  xc  \2}  =  p  +  cj2  (Bll) 
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To  find  the  variance  of  zx ,  we  first  find, 

E{zl)  =  E{ |  x,  +  x,  |4}  =  E{ |  x,  |4}  +  E{ |  x,  |4}  +  4E{\ x,  |2}£{|  xc  |2}  (B12) 

E{\  xt  |4}  =  2er4  (B13) 


E{  |xj4} 


2(v  +  l)  2 

fl 

V 


(B14) 


Inserting  (B11)-(B14)  into  (BIO),  we  have, 

U  =  v(l+/?)2  (B15) 

where  /?  =  <r2  ///  is  the  SCR  of  the  data.  Equation  (B15)  is  identical  to  what  is  given  by 
Watts  (1987)4. 


Similarly,  the  third  normalised  moment  of  K-distribution  is, 
6(v1+2)(v1+l)_  E{z\) 

vi  E3{Zl} 

After  some  elementary  algebra,  (B16)  becomes, 

3vj  +  2  _  (v  +  2)(v  + 1)  /  v2  +  3 /?(v  + 1)  /  v  +  3/?2  +  ft 

v,2  (1  +  ^)3 


(B16) 


(B17) 


When  ^  0 ,  (B17)  is  a  quadratic  equation  of  the  unknown  y  which  is  solvable  once  the 
shape  parameter  of  clutter,  v ,  and  the  SCR,  j3 ,  are  given.  Its  explicit  solution  is. 


vi 


3  +  V9  +  8  A 
2  A 


(B18) 


where  A 


Ml  +  P)  +  2 

vQ-  +  Pf 


We  thus  have  derived  two  Formulae,  i.e.,  (B15)  and  (B18)  for  estimating  the  shape 
parameter  v .  If  the  distribution  of  target  embedded  in  clutter  is  exactly  a  K-distribution, 
the  two  formulae  should  result  in  an  identical  value.  Since  the  true  distribution  of  the 
combined  two  is  unknown,  and  the  K-distribution  is  only  an  approximation,  the  shape 


4  It  is  noted  that  in  the  Watts'  paper  (1987),  the  clutter-to-noise  ratio  (CNR)  is  used  whereas  here  the 
signal-to-clutter  ratio  (SCR)  is  used. 
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parameter  estimated  by  (B15)  and  (B18)  may  be  slightly  different.  Often(B18)  results  in  a 
slightly  smaller  value,  and  the  corresponding  distribution  may  fit  marginally  better  in  the 
tail  region. 

The  accuracy  of  using  a  new  K-distribution  to  approximate  the  distribution  of  Gaussian 
target  signals  embedded  in  K-distributed  clutter  may  be  examined  by  comparison  between 
the  pdf  numerically  integrated  by  (B9)  (treated  as  the  'true'  distribution  when  the 
integration  interval  is  sufficiently  small)  and  the  pdf  of  the  K-distribution  whose  shape 
parameter  is  determined  by  (B15)  or  (B18).  Figure  17  shows  such  a  comparison.  It  can  be 
seen  from  the  figure  that  the  K  approximation  is  almost  identical  to  the  true  distribution, 
though  the  upper  tails  of  the  two  separate  apart  slightly  at  the  far  end.  However,  the 
accuracy  of  using  a  new  K  distribution  to  approximate  the  distribution  of  Gaussian  target 
signals  embedded  in  K-distribution  clutter  would  be  sufficient  for  the  purpose  of  radar 
detection  analysis. 


Intensity  (dB) 

(a)  pdf  on  linear  scale 


(b)  pdf  on  log  scale 


Figure  17:  Pdf  comparison  between  numerically  integrated  (viewed  as  'true  distribution')  and  the  K 
approximation  (v  =  l.O,  /u  =  1.0  and  SCR  =  5dB  )  The  shape  parameter  v1  is 
determined  by  (B18). 


The  Monte  Carlo  simulation  was  further  used  to  examine  the  agreement  between  the 
assumed  distribution  and  the  actual  distribution  of  target  plus  clutter.  In  the  simulation,  10 
million  ( 107 )  samples  were  generated,  and  the  histograms  of  data  were  compared  with  the 
pdf  of  the  K-distribution.  Figure  18  and  Figure  19  show  the  goodness-of-fit  of  the 
histograms  of  data  and  the  pdfs  of  the  corresponding  K-distributions.  It  can  be  seen  that 
agreement  between  the  two  is  superior,  and  hence  the  distribution  of  a  Gaussian  target 
embedded  in  K-distributed  clutter  can  be  confidently  modelled  as  a  new  K-distribution. 
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Intensity  (dB) 

(b)  Pdf  on  log  scale  to  show  the  goodness-of-fit  for  tails 


Figure  18:  Goodness-of-fit  for  distribution  of  Gaussian  target  embedded  in  K-distributed  clutter 
(ju  =  1.0,  v  =  1.2,  SCR  =  3dB). 
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Intensity  (dB) 

(b)  Pdf  on  log  scale  to  show  the  goodness-of-fit  for  tails 

Figure  19:  Goodness-of-fit  for  distribution  of  Gaussian  target  embedded  in  K-distributed  clutter 
(ju  =  1.0,  v  =  1.2,  SCR  =  \0dB ). 


It  is  found  that  once  a  Gaussian  target  is  embedded  in  the  K-distributed  clutter,  the 
distribution  of  the  combined  two  becomes  much  less  spikier  than  the  original 
K-distributed  clutter,  i.e.,  the  new  shape  parameter  vl  is  usually  larger  than  v .  If  SCR  is 
large,  the  combined  distribution  approaches  Gaussian  (the  intensity  approaches 
exponential),  as  the  combined  signal  has  a  dominant  Gaussian  component.  For  numerical 
calculations,  the  K-distribution  with  a  shape  parameter  >  50  can  be  approximated  as  a 
Gaussian  distribution  (exponential  distribution  for  the  intensity).  Figure  20  shows  the 
shape  parameter  vt  as  a  function  of  v  and  SCR.  If  can  be  see  that  for  moderate  sea  clutter 
having  a  shape  parameter  of  5  or  larger,  the  distribution  of  Gaussian  target  plus  clutter  can 
be  treated  as  Gaussian  distributed  once  the  SCR  is  greater  than  3dB.  However,  for  spikier 
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sea  clutter  with  a  shape  parameter  less  than  1,  the  distribution  of  the  combined  two  slowly 
approaches  Gaussian  only  when  SCR  close  or  greater  than  lOdB. 


Figure  20:  Shape  parameter  vl  rapidly  increases  with  the  increase  in  SCR  (solid-line  and  broken- 
line  curves  are  calculated  by  (B15)  and  (B17),  respectively). 


Another  way  to  examine  whether  a  distribution  is  exponential  is  to  examine  its  statistical 
properties  including  mean,  median,  variance,  skewness  and  kurtosis.  Table  A1  lists  these 
statistical  measures  for  the  Monte  Carlo  simulated  10  million  (107)  samples  of  Gaussian 
target  embedded  in  the  K-distributed  clutter.  It  shows  that  with  an  increase  in  SCR,  the 
distribution  of  the  combined  signal  intensity  indeed  approaches  exponential.  Figure  21 
examines  the  goodness-of-fit  for  a  case  of  Gaussian  target  embedded  in  the  K-distributed 
clutter,  with  parameters  of  p  =  1.0,  v  =  5.0,  SCR  =  4dB .  It  shows  that  the  distribution  is 
almost  identical  to  the  exponential  distribution. 


Table  Al:  Statistical  measures  of  Monte  Carlo  simulation  (10b  samples)  in  comparison  with 
theoretical  values  of  an  exponential  distribution. 


Mean 

Median 

Variance 

Skewness 

Kurtosis 

Exponential  distribution 
p(z)  =  exp  (-z/  jU)/  p 

M 

(1.000) 

//  ln(2) 
(0.6931) 

2 

jU 

(1.000) 

2 

6 

10  million  samples  of 
combined  target  +  clutter, 

v  =  5,  SCR  =  4 dB 

1.000 

0.686 

1.032 

2.097 

6.821 

10  million  samples  of 
combined  target  +  clutter, 

v  =5,  SCR  =  lOdB 

1.000 

0.692 

1.004 

2.014 

6.107 
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Intensity  (dB) 

(b)  Pdf  on  log  scale  to  show  the  goodness-of-fit  for  tails 

Figure  21:  Goodness-of-fit  for  distribution  of  Gaussian  target  embedded  in  K-distributed  clutter 
( ju  =  1.0,  v  =  5.0,  SCR  =  AdB ).  The  approximated  K-distribution  is  almost  identical 
to  the  exponential  distribution. 
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Appendix  C:  Temporal  and  Spatial  Correlations  of 
Compound  K-Distributed  Clutter 

The  compound  K-distribution  assumes  that  the  data  consists  of  a  fast-varying  component 
modulated  by  a  slowly-varying  component.  The  fast-varying  component  commonly  refers 
to  the  speckle  component  that  is  a  Gaussian  random  process  with  zero  mean  and  unit 
variance.  The  slowly-varying  component  is  also  called  underlying  mean  whose  intensity 
has  a  gamma  distribution.  Because  of  the  nature  of  these  two  components,  caution  has  to 
be  exercised  when  estimating  temporal  correlation  and  spatial  correlation. 

C.l.  Temporal  Correlation  (Correlation  in  Azimuth) 

It  often  needs  to  know  the  temporal  correlation  between  data  samples  collected  by  a  pulse 
train  in  a  coherent  processing  interval  (CPI).  Often  we  can  assume  during  the  CPI,  the 
slowly-varying  component  remains  unchanged,  i.e.,  the  slowly  varying  component  is  fully 
correlated  (this  is  how  the  compound  K-distributed  clutter  got  its  name)5.  The  covariance 
matrix  of  the  data  is  written  as, 

M=^{xx//}  (Cl) 

where  expectation  is  with  respect  to  the  pulse  for  the  temporal  correlation,  and 
x  =  [x[0]  •  •  •  x[N  -  l]f  is  a  N  xl  vector  collected  by  N  pulses.  Each  measurement  may 

be  further  written  as  a  product  of  fast-varying  and  slowly-varying  components,  according 
to  the  model  assumption,  as, 

x[n ]  =  yfr xf[n ]  (C2) 

where  xf[n ]  is  the  fast-varying  component  which  is  a  complex  Gaussian  variable,  r  is  the 

underlying  mean  which  is  constant  in  a  CPI.  Because  two  components  are  independent, 
we  have, 

M  =  E\xxH}=E{r}E^ifxHf}= /uMf  (C3) 

where  M;  =  e\x  /  x///  }  is  the  covariance  matrix  of  the  fast-varying  component  and  ju  the 

clutter  mean.  Therefore,  the  temporal  correlation  can  be  estimated  using  data  samples  and 
the  correlation  of  fast-varying  component  is  just  the  covariance  matrix  normalised  by  its 
mean.  Under  the  assumption  of  wise  sense  stationary,  M  f  has  a  Toeplitz  structure  of. 


5  With  reference  to  a  maritime  search  radar,  there  is  high  superposition  of  antenna  footprints  with 
respect  to  successive  pulses  in  a  CPI,  so  the  texture  of  sea  clutter  can  be  assumed  to  be  completely 
correlated. 
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A  ‘  ‘ '  Av-i 
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i  •••  ; 

:  a 

* 

Av-i 

»  , 

•••  A  1 

where  pk  =  —  is{x[7]x*[z  +  k]}r 


i,  k  =  0, 1,  •  •  • ,  N  —  1 . 


(C4) 


C.2.  Spatial  Correlation  (Correlation  in  Range) 

The  correlation  of  texture  may  be  found  from  the  spatial  correlation.  However  in  order  to 
find  the  correlation  of  texture,  we  normally  need  to  manipulate  data  in  the  intensity  (or 
amplitude)  domain  rather  than  the  complex  (IQ)  domain.  Because  if  we  consider  the 
correlation  between  range  bin  i  and  range  bin  i  +  k ,  k  ^  0 6 , 

E\^T[i]xf[i]^T[i  +  k]x*f[i  +  &]}  =  e\J  r[i]r[i  +  k]^E  {x  f[i]x*f[i  +  /r]}  (C5) 

However,  if  the  interval  between  bin  i  and  bin  i  +  k  is  greater  than  the  radar  range 
resolution,  e\x f[i\x*f[i  +  &]}=  0 .  Therefore,  even  if  E |yV|7]r[/'  +  /Q } ^  0 ,  its  value  is  not 

measurable  using  the  IQ  data.  Therefore,  the  correlation  of  texture  has  to  be  found  using 
the  intensity  (or  amplitude)  data.  First  we  define  the  correlation  coefficient  of  texture,  p , 
and  the  correlation  coefficient  of  intensity,  rj ,  respectively,  as. 


E{r[i]r[i  +  k\}-E2{r} 

Pk  ~  ,  .  K  ~  A  b  ’  ’  ’ 

var(r) 

(C6) 

_E{z[i]z[i  +  k]}-E2{z} 

Ik  ~  ,  .  A,  -  U,  1,  •  •  • 

var(z) 

(C7) 

Obviously, 

rj0  =  Po  =  1 

(C8) 

Below  we  will  find  the  relationship  between  rj k  and  pk  for  k  ^  0  . 
Denoting  z  =\x\2  =  t\  x  f  2  =  r  z,,  we  have. 


6  In  order  to  simplify  symbols,  the  same  index  notation  is  used  to  represent  either  temporal  series 
(the  index  refers  to  pulse  numbers)  or  spatial  series  (the  index  refers  to  range  bin  number).  There 
should  be  no  confusion  under  the  context. 
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E{z[i]z[i  +  A']}  =  £'{r[/]r[/  +  k]zf\i]zf\i  +  A:]} 


(C9) 


Since  the  fast-varying  and  slowly  varying  components  are  independent,  we  have, 

£'{r[/]r[/  +  k]zf[i]zf[i  +  A]}=  As{r[/]r[7  +  k]}E\zf[i]zf[i  +  A]}  (CIO) 


For  A  ^  0 ,  the  value  of  E \z  f  [ i]z  f  [l  +  A] }  can  be  calculated  by  the  Isserlis'  Theorem  (see 

Michalowicz,  J.V.,  Nichols,  J.M.,  Bucholtz,  F.,  and  Olson,  C.C.,  "An  Isserlis'  theorem  for 
mixed  Gaussian  variables:  application  to  the  auto-bispectral  density",  ].  Statistical  Phys., 
2009, 136,  (1),  pp.  89-102,  for  details),  as. 


E{zf[i]zf[i  +  A]}  =  E^cf[i]x*  \}\xf\i  +  k]x*f[i  +  A]] 

=  1  +e\x f\j\x f\i  +  k~^E\xf[i\x*f\i  +  k~^+E\x f[i\x*f\i  +  k~^E\xf[i\x f[i  +  k~^ 


for  k  ^  0 


(Cll) 


If  the  interval  between  bin  i  and  bin  i  +  k  is  greater  than  the  radar  range  resolution,  the 
last  two  items  of  (Cll)  become  zero,  so  E \z f[i]z f[i  +  k]\  =  1,  the  above  correlation 
simplifies  to, 

E{c[i\u[i  +  A]}  =  E{z[i]z[i  +  A]}  for  A  ^  0  (C12) 


Combining  (C6),  (C7)  and  (C12),  and  noticing  E{t)  =  E{z]  =  //,  one  achieves. 


nk 


var(r) 

var(z) 


A  =1,2--- 


(C13) 


Equations  (C8)  and  (C13)  indicate  once  we  know  the  correlation  of  z  or  r,  the  correlation 
of  the  other  can  be  determined. 


The  texture  r  is  gamma  distributed,  and  the  intensity  z  is  single-look  or  multi-look 
(multi-pulses)  K-distributed,  their  variances,  respectively,  are. 


var(r)  =  /u2  !  v 


var(z)  = 


V  +  \  +1  2 


Nv 


A 


(C14) 

(CIS) 


where  N  is  the  number  of  multi-looks  (multi-pulses).  Therefore,  because  of  the  effect  of 
the  fast-varying  component  that  is  uncorrelated  and  randomly  varies,  the  correlation  of 
the  intensity  z  is  generally  weaker  than  the  correlation  of  the  texture  r  that  is  the 
originator  of  the  correlation.  Only  if  the  number  of  multi-looks  reaches  infinity  when  the 
fluctuations  of  the  fast-varying  component  disappears  (averaged  to  its  mean  value  for 
every  range  bin),  two  correlations  become  identical. 
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Inserting  (C14)  and  (C15)  into  (03),  we  have. 


vk  =  Pk 


N 

v  +  N  +  l 


k=  1,2- 


(06) 


The  correctness  of  the  above  derivation  is  confirmed  by  the  Monte  Carlo  simulation.  In  the 
simulation,  we  first  generated  the  slowly  varying  component,  gamma,  that  has  a  shape 
parameter  of  v  =  1.2  and  a  correlation  coefficient  of, 

pk  =[0.7  +  0.3cos(0.12;r/t)]e-A/12  k  =  0,1,  -  (07) 


The  multi-look  K  distributed  data  was  then  generated  by  modulating  the  uncorrelated 
multi-look  fast-varying  Gaussian  component  to  the  correlated  gamma.  The  correlation 
coefficient  for  the  simulated  multi-look  K  data  was  regressed  and  compared  to  the 
theoretical  value,  and  details  are  shown  in  Figure  22.  It  can  be  seen  that  the  simulated 
results  match  the  theory.  Only  when  the  number  of  multi-looks  in  the  averaging 
processing  (multi-pulse  averaging  processing)  becomes  large,  the  correlation  of  the  multi¬ 
look  data  approached  the  correlation  of  the  slowly-varying  component. 


(a)  Correlation  of  gamma  (b)  correlation  of  K 

Figure  22:  (a)  the  designated  correlation  of  the  gamma  component  and  (b)  the  correlation  of  the 
multi-look  K  data  (the  cases  of  1-look  and  16-looks  are  shown). 


For  range  over-sampled  data,  when  the  range  interval  between  bin  i  and  bin  i  +  k(k*  0) 
is  smaller  than  the  range  resolution,  the  fast-varying  component  is  also  correlated  in  range, 
resulting  in  E\zf[i]zf[i  +  &]}>  1  (see  (Cll)).  So  the  correlation  of  the  intensity  z  will  be 
jointly  contributed  by  the  fast-varying  component  and  the  slowly-varying  component. 
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Appendix  D:  Inverse  of  Covariance  Matrix 


Suppose  a  covariance  matrix  is  symmetrical,  positive  definite  and  has  a  Toeplitz  structure 
of. 


1 

A  Pi 

Pn- i 

* 

A 

1  Px  ■■■ 

Av-2 

* 

Pn- 2 

•••  A  1 

A 

* 

Pn- i 

*  * 

•••  Pi  A 

1 

(Dl) 


The  inverse,  while  still  be  symmetrical  and  positive  definite,  does  not  in  general  have  a 
specific  structure. 


However,  if  pk  k  =  0,  •  •  • ,  N  - 1  is  a  geometric  series,  e.g.. 


M  = 


1 

* 

P 


( PN~2r 

{pN-ir 


i  p 

*  i 

••  p  1 

••  iP  )  P 


P 

P 


N- 1 

N-2 


P 

1 


(D2) 


Then  the  inverse  of  M  has  a  simple  structure  as. 


M_1 


a  c  0  0 

c  b  c  0 

0  c  b  c 
0  •••  0  c  a 


where 


a  = 

<b  = 


1 


1-1  P 
1+1  P 

1-1  P 


|2 

2 

|2 


C  = 


~P 


1-1  P 


2 


(D3) 


(D4) 
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T o  prove  (D3),  we  denote  Al  .  for  element  of  A ,  i,  j  =  1,  •  •  •  N ,  and 


1 

P  P2 

...  pN~ 1 

c  0 

0" 

* 

p 

1  P 

N—2 

...  p 

* 

c 

b  c 

0 

(P  ) 

* 

...  p 

1  P 

0 

•••  c*  b 

c 

1 

* 

-  (P2)* 

P  1  J 

0 

■■■  0  c* 

a 

It  is  not  difficult  to  write, 

Au  =a  +  pc*  =1,  ANN  =  p*  c  +  a  =  1 
An=  p*c  +  b  + pc  =  1,  i  =  2,---,N-l 

Aj  j  =  pj~'~lc  +  p’-'b  +  pH+xc  =  pj~'~x  (c  +  pb  +  p~~ c *  )=  0,  i  <  j 
A.  .=  A*  .=  0  i  >  j  (A  is  Hermittian) 


(D5) 


(D6) 
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Appendix  E:  Distribution  of  Muli-look  Correlated  K- 

Distributed  Data 


Multi-look  processing,  coherently  or  non-coherently,  is  a  common  practice  for  radar  data 
processing.  In  some  cases,  data  are  correlated.  We  want  to  derive  the  distribution  (pdf  and 
cdf)  of  multi-look  correlated  data,  including  Gaussian  and  K-distributed  data. 

Denoting  M  the  normalised  covariance  matrix  of  correlated  complex  Gaussian  vector  x, 
the  correlation  can  be  expressed  by  spherically  invariant  random  process  (SIRP) 
(Rangaswamy  et  al.  1993;  Antipov  1998).  Letting  t  be  a  mutually  uncorrelated  circularly 
complex  Gaussian  vector  with  vector  zero  mean  and  covariance  matrix  of  c2I  (I  is  an 
elementary  matrix),  the  correlated  Gaussian  vector  x  may  be  written  as, 

x  =  M1/2t  (El) 

where  M1/2  =  EA1/2E//,  E  and  A  are  eigenvectors  and  eigenvalues  of  M  (note  that  M 
is  a  Hermitian  matrix,  e.g.,  symmetrically  positive  definite  matrix,  and  in  addition  it  has  a 
Toeplitz  structure).  Obviously  the  covariance  matrix  of  x  is, 

£{xx"  }=  £{M1/2tt"  (M1/2)"  }=  a2M  (E2) 


So  the  normalised  covariance  matrix  is  M  =  — E |xxw ] ,  and  =  E{\  x\2}  .  Now  we 

<7 

1  2  1 

want  to  derive  the  pdf  of  — x\i\  =  —  x  x  .  First  we  have, 
r  N  N 


—  x"x  =  —  (M1/2t)"  (Ml/2t)  = 
N  N 


1  N 

— Z4K 

nU 


N 

=  Wtt 


(E3) 


where  kt  =Ai/N  and  Ajr  i  =  1 , •  •  • ,  /V ,  are  the  eigenvalues  of  M.  Note  that  i  =  1 , •  •  • ,  /V , 
are  mutually  uncorrelated  and  each  is  a  complex  Gaussian  distributed,  and  ]  tt  |2  is 
exponentially  distributed  and  has  a  mean  of  cr2 . 

1  2 

According  to  (E3),  the  pdf  of  —  I  Aj\  T  is  identical  to  the  pdf  of  X/=i  \  t[i]  |2 .  Since  tt 

( i  =1,  •  •  N )  are  mutually  uncorrelated,  the  pdf  of  X/=i k,  |  /[/]  |2  can  be  derived  from  the 
so-called  convolution  formula  (Wilks,  1962,  p.  204)  of. 


PN(t)  =  jp(t  -x)pN_!  (x)dx 
0 


(E4) 
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where  p(x)  =  px  (x) .  Using  (E4),  we  derive  the  pdf  of  t 


N 


N ' 


as. 


Pn(  0 


1  ^k^exp^t/jk^2)) 

2  2-u  N 

^  r Wt-K) 

m=  1 
m^i 


The  cdf  is. 


N  pt 

Ar(0  =  !-Z  — 


_1  exp(-f /(£;cr)) 

Wk,  -K) 


The  nth  moment  of  t  is. 


k 


N+n- 1 


i m-kj 


m= 1 
m&i 


(E5) 


(E6) 


(E7) 


It  is  known  that  the  multi-look  uncorrelated  Gaussian  data  has  a  gamma  distribution  of. 


p{z)  = 


r(A)((j2  in) 


j*  zN~l  exP| 


V  a  I N  j 


(E8) 


1  Jv 

where  z  =  — x[n]  \  ,  and  N  is  the  number  of  multi-looks. 

N  ;=i 


The  distribution  of  multi-look  correlated  Gaussian  data  is  given  by  (E5)  which  is  not  a 
gamma  distribution.  However,  for  simplification,  it  may  be  approximated  by  a  gamma 
distribution  whose  mean  and  equivalent  number  of  independent  looks  may  be  found  by 
equating  the  mean  and  the  variance  of  the  both  distributions.  In  general,  the  equivalent 
number  of  independent  looks  is  greater  than  1  (fully  correlated  case)  and  less  than  N 
(fully  uncorrelated  case).  For  example,  if  the  covariance  matrix  has  a  Toeplitz  structure,  as. 


1 

A  Pi 

Am 

* 

A 

1  a  ••• 

Pn-1 

* 

Pn- 2 

•••  A  1 

Pi 

* 

Pn-  i 

*  * 

Pi  A 
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(E9) 
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By  equating  the  first  and  second  moments  of  (E5)  and  (E8),  we  found  the  equivalent 
number  of  independent  looks  to  be, 

k  = - - - - -  N>  2  (E10) 

2  N-l  v  7 

l  +  -X(N-n)\pn  |2 
N  n=  1 

Figure  23  compares  the  exact  distribution  and  the  approximated  gamma  distribution  for 
multi-look  of  Gaussian  distributed  data  with  N  =  8  ,  and  correlation  coefficients  of 
pn  =  exp(-n) .  It  can  be  seen  that  that  approximated  gamma  distribution  starts  deviating 
from  the  exact  distribution  in  the  far  upper  tail  region. 


(a)  pdf  ( pn  =  exp(-«) )  (b)  1-cdf  ( pn  =  exp(-») ) 

Figure  23:  Distribution  of  t  =  AG'ZmI  x\i\  f  is  approximated  as  a  gamma  distribution  when 
samples  are  correlated.  The  blue  dots  represent  data's  distribution  and  the  green  lines 
stand  for  the  approximated  gamma  distribution  (N  =  8,  pn  =  exp(-«)  which  result  in 
k  =  6.3102). 


The  distribution  of  multi-look  correlated  K-distributed  data  can  be  derived  in  the  similar 
way.  We  assume  that  the  slowly-varying  component  remains  constant  during  the  multi¬ 
look  process,  i.e.,  the  slowly-varying  component  is  fully  correlated  during  the  period  of 
multi-looks  and  the  fast-varying  component  is  partially  correlated  and  has  the  covariance 
matrix  of  M .  The  corresponding  pdf  and  cdf  are  given  by  (Ell)  and  (E12),  respectively. 
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where  z  =  N  'Xilil  Af\  |2  .  Similarly  the  exact  distribution  given  by  (Ell)  can  be 
approximated  by  a  multi-look  K-distribution  given  by  (B6)  with  an  equivalent  number  of 
independent  looks.  Figure  24  compares  the  exact  distribution  and  the  approximated 
distribution  for  multi-look  of  K-distributed  data  with  N  =  8 ,  and  correlation  coefficients 
of  pn  =  exp(-n)  .  Although  the  correlation  is  the  same  as  the  above  example,  the 
approximated  pdf  now  demonstrates  much  better  agreement  with  the  exact  pdf.  Because 
for  the  K-distribution,  its  upper  tail  is  primarily  determined  by  the  shape  parameter  (the 
slow-component),  and  the  approximation  processing  only  affects  the  fast-varying 
component,  but  has  no  effect  on  the  slowly-varying  component.  As  a  result,  the 
approximation  for  the  K-distribution  is  better  for  small  shape  parameters  than  larger  ones. 
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Normalised  intensity  (dB) 


(b)  pdf  on  log  scale 


Normalised  intensity  (dB) 

(c)  1-cdf  on  log  scale 

Figure  24:  An  example  of  multi-correlated-look  K-distribution,  N  =  8,  pn  =  e~"  and  v  =  \.2.  The 

equivalent  uncorrelated  look  number  is  *:  =  6.3102  for  the  approximate  multi-look  K 
distribution. 
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